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INTRODUCTION 


It  is  inevitable,  however  unfortunate,  that  any  large 
collection  of  anthropometric  data  will  contain  an  unknown  number 
of  errors.  While  the  sources  of  such  mistakes  are  many,  their 
chief  cause  is  human  fallibility.  In  the  1967  U.  S.  Air  Force 
Anthropometric  Survey,  for  example,  2420  flight  crewmen  were 
measured  for  190  variables — an  awesome  total  of  459,800  meas- 
urements. It  requires  but  a momentary  mental  distraction 
during  any  of  the  measuring  or  data  processing  steps  to  cause 
an  error  to  slip  into  the  recorded  material.  Mistakes  can 
occur  during  the  marking  or  positioning  of  subjects,  during  the 
measuring  of  subjects  or  the  reading  of  instruments,  as  well  as 
during  the  recording,  transposition  and  punching  of  data. 

It  is  of  some  importance,  therefore,  to  introduce  methods 
for  the  checking  and  evaluation  of  anthropometric  data  which 
can  provide  some  degree  of  assurance  that  errors  of  measurement 
and  recording  are  identified  and  reduced  to  a minimum.  This 
is  accomplished  at  increasingly  more  sophisticated  levels  with 
the  aid  of  a high-speed  computer  which  can  rapidly  check  and 
evaluate  vast  numbers  of  numerical  observations  in  a systematic 
fashion . 

The  two  editing  programs  described  here  were  developed  by 
Edmund  Churchill  of  the  Anthropology  Research  Project  under 
contract  to  the  Aerospace  Medical  Research  Laboratory  and,  when 
used  together,  effect  a "coarse"  and  "fine"  sifting  of  the  data. 
The  XVAL  (extreme  VALues)  routine  scans  for  gross  errors  by 
sorting  out  the  ten  highest  and  ten  lowest  values  from  a given 
group  for  examination.  The  EDIT  (EDITing)  program  is  used  to 
test  each  data  point  for  each  subject  by  comparing  the  measured 
value  with  a predicted  value  obtained  from  regression  equations 
and  flagging  those  which  exceed  set  limits.  While  either  one 
of  the  two  programs  can  be  used  without  the  other,  they  have 
been  developed  as  complementary  procedures  and  are,  by  far,  more 
effective  together  than  they  are  separately. 

In  general,  the  editing  of  anthropometric  data  consists  of 
three  sequential  steps: 

• selection  of  a suitable  test  function, 

• selection  of  an  appropriate  acceptance/rejection  level, 

• selection  of  a course  of  action  when  rejection  of  a data 
point  occurs. 

Both  editing  programs  described  here  routinely  achieve  the 
first  objective,  leaving  to  the  user  (in  the  case  of  the  EDIT 
program)  only  the  task  of  selecting  the  variables  to  be  used  in 
the  regression  equations. 
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Selection  of  an  appropriate  acceptance/rejection  level 
(second  step)  is,  of  course,  left  to  the  user  and  can  be  as 
stringent  or  as  lax  as  desired.  The  user  must  decide  on  the 
trade-offs  between  establishing  a low  rejection  level  resulting 
in  the  increased  possibility  of  retaining  erroneous  values,  and 
setting  very  stringent  limits  thereby  increasing  the  possibility 
of  rejecting  valid  data  because  they  deviate  in  an  insignificant 
way  from  the  rest  of  the  distribution  of  values. 

The  third  step  is  also  left  to  the  user's  discretion.  When 
a data  point  is  obviously  in  error,  the  decision  to  correct, 
delete,  substitute  or  remeasure  is  one  which  must  be  made  by 
the  user  based  on  his  knowledge  and  judgment  concerning  the 
worth  of  the  data  point.  If  the  faulty  value  is  one  of  a large 
series  and  does  not  affect  the  sample  distribution,  deletion  or 
substitution  may  be  in  order.  If  the  data  point  is  important 
to  the  distribution,  then  remeasurement  may  be  called  for  if 
an  appropriate  substitution  cannot  be  found. 

While  the  techniques  discussed  here  and  more  fully  described 
on  the  ensuing  pages  are  designed  and  most  often  used  for  editing 
anthropometric  survey  data  when  both  the  sample  size  and  the 
number  of  variables  are  relatively  large  (n  >100) , they  can  be 
applied  equally  well  to  studies  involving  fewer  subjects  and 
variables . 
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THE  EXTREME  VALUES  (XVAL)  PROGRAM* 


The  basic  purposes  of  the  extreme  VALues  (XVAL)  program  are 
to  provide  an  initial  scanning  of  a set  of  data,  singling  out 
extreme  outlyers  worthy  of  examination  for  errors  in  recording, 
punching,  or  coding,  and  to  provide  first  approximations  for  the 
means  and  standard  deviations  for  each  variable  in  the  total 
sample.  This  program  is  the  first  one  used  in  editing  our 
anthropometric  data.  We  check  the  outlyers,  correcting  or 
discarding  those  in  obvious  error.  Usually  the  data  are  then 
rerun  and  the  process  repeated  until  the  values  at  the  tails  of 
the  distributions  are  deemed  acceptable.  For  many  sets  of  data 
this  final  output  may  provide  all  the  editing  needed. 

The  XVAL  program  evolved  from  a series  of  rather  simple 
programs  designed  to  provide  a listing  of  the  very  smallest 
and  the  very  largest  values  for  each  of  several  sets  of  data. 

As  the  program  has  evolved,  the  present  output  also  provides 
many  of  the  standard  summary  statistics,  two  optional  listings 
and  an  optional  set  of  punched  cards,  as  well  as  a number  of 
other  pieces  of  information  all  of  which  we  shall  describe. 


THE  OUTPUT 

The  primary  XVAL  output  is  shown  in  Figure  1.  The  values  in 
each  column  refer  to  the  variable  whose  name  and  number  are 
listed  at  the  top  of  the  column.  Thus,  the  first  column  contains: 

(1)  The  ten  smallest  values  of  age. 

(2)  Subject  or  record  number.  (As  can  be  seen,  subject 
#120  had  a recorded  age  of  13.5  years,  clearly  a 
suspect  value  in  a military  survey.) 

(3)  The  ten  largest  values  of  age. 

(4)  The  arithmetic  mean  of  the  entire  sample. 

(5)  The  standard  deviation  of  the  entire  sample. 

(6)  The  coefficient  of  variation  for  the  entire  sample. 

(7)  A value  expressing  the  variation  within  the  10  top 
values  as  a fraction  of  the  difference  between  the 
10th  smallest  and  the  10th  largest  values. 


* A complete  printout  of  the  XVAL  program  is  included  as  Appendix  A. 
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Figure  1.  The  primary  output  of  the  XVAL  program  (age  in  tenths 

of  years,  weight  and  grip  strength  in  tenths  of  kilograms 
values  for  all  other  variables  in  millimeters) - 


(8)  A value  expressing  the  variation  within  the  10  bottom 
values  as  a fraction  of  the  difference  between  the 
10th  smallest  and  the  10th  largest  values. 

(9)  Symmetry  or  skewedness  (veta  I) . 

(10)  Kurtosis  or  peakedness  (veta  II) . 

(11)  Estimate  of  the  mean  after  elimination  of  the  10  top 
and  10  bottom  values. 

(12)  Estimate  of  the  standard  deviation  after  elimination 
of  the  10  top  and  10  bottom  values. 

(13)  The  difference  between  the  mean  (4)  and  the  estimated 
mean  (11)  expressed  as  a percent  of  the  estimated 
standard  deviation. 

(14)  The  difference  between  the  standard  deviation  (5)  and 
the  estimated  standard  deviation  (12)  expressed  as 

a percent  of  the  estimated  standard  deviation. 

(15)  The  number  of  data  points  for  this  variable. 

Note  that  if  the  sample  size  for  any  variable  is  less  than 
30,  the  n-20  mean  and  standard  deviation  estimates  and  the  sta- 
tistics TOP  and  BOT  are  not  computed  for  that  variable.  If 
the  sample  size  is  less  than  or  equal  to  20,  only  the  total 
sample  mean  and  standard  deviation  are  computed  along  with 
some  range  card  values. 

The  label  (16)  at  the  bottom  of  the  page  in  the  center,  and 
the  date  (17)  at  left  have  been  read  in  from  cards.  The  page 
number  (18)  at  right  is  kept  track  of  by  the  program. 

A summary  table  (Figure  2)  follows  this  output.  It  contains 
for  each  variable,  (1)  the  mean  value,  (2)  the  standard  deviation, 
(3)  and  (4)  the  measures  of  skewness  and  kurtosis,  (5)  the  coeffi- 
cient of  variation,  (6)  and  (7)  the  percentage  difference  for  the 
mean  and  standard  deviation  as  described  in  13  and  14  above,  (8) 
the  sample  size  for  that  variable,  and  (9)  the  minimum  value  for 
that  variable.  The  final  seven  values  contained  in  the  summary 
statistics  for  each  variable  make  up  the  optional  card  output 
(Figure  3)  for  XVAL  along  with  the  variable  number  and  name 
as  follows: 

The  Summary  Statistics  The  Optional  Punched  Name-Range  Cards 

Columns  1-4:  the  variable  number 

Columns  7-24:  the  variable  name 
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Figure  3. 


The  optional  name-range  card  output 
for  the  XVAL  program. 
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The  Summary  Statistics  The  Optional  Punched  Name-Range  Cards 
(continued) 

(10)  MIN  - ' Columns  25-32:  a value  generally 

slightly  smaller  than 
the  minimum 


(11)  MAX  - 

(12)  AVG  - 


Columns  33-40:  the  maximum  value 

Columns  41-48:  an  approximate  mean 

value 


(13)  INTVl 


Columns  49-54:  a suggested  interval 

width  for  frequency 
tables  which  are  lim- 
ited to  a maximum  of 
30  intervals 


(14)  INTV2  - Columns  55-60:  a suggested  interval 

width  for  frequency 
tables  which  are  lim- 
ited to  a maximum  of 
50  intervals 


(15)  CF-1 


Columns  61-70:  a constant  for  con- 
verting the  data  for 
that  variable  to  the 
usual  metric  units 


(16)  CF2  - Columns  71-80:  a constant  for  con- 

verting the  data  for 
that  variable  from  the 
usual  metric  units  to 
the  usual  English 
units 


At  the  bottom  of  each  page  of  summary  statistics,  we  have 
again  the  date,  the  label  for  this  run,  and  the  page  number. 

Besides  the  optional  card  output,  there  is  also  optional 
printout  material  consisting  of  a listing  of  missing  data  values; 
i.e.,  values  read  in  as  zero  (Figure  4a),  and  a listing  of  miss- 
ing subject  (record)  numbers  (Figure  4b) . The  latter  actually 
results  from  each  instance  in  which  a subject  number  read  in  is 
not  precisely  one  more  than  the  previous  subject  number.  Data 
records  not  in  order,  or  ones  with  duplicate  subject  numbers  will 
trigger  listings. 

The  only  other  output  of  program  XVAL  itself  is  printing  of 
the  namelist  XVAL’s  variable  values  and  messages  associated  with 
these  values  (Figure  5) . These  values  are  discussed  in  detail 
beginning  on  page  26. 
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NR 
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NQ 

= 25, 

XKST 

= • 1E+00 , 

YKST 

= • 3 9370  0 SE  + Q 0 , 

IP 

II 
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NUNIT 

•s 

J- 

II 

SENO 
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■-MISSING  VALUE  LIST  WILL  BE  PRINTED 

MS.GT.O- 

■-NON -SEQUENTIAL  SUBJECTS  WILL  BE  LISTED 

IP. LE. 4— FRACTIONAL  INTERVALS  POSSIBLE 

NR.GT.O— NAME  RANGE  CARDS  WILL  BE  PUNCHED 
MAXIMUM  NUMBERS  OF  INTERVALS  WILL  BE,  45  £ 25 

Figure  5.  The  namelist  XVAL. 


15 


Other  outputs  attributable  to  subroutines  used  with  XVAL  are, 
in  order  of  appearance,  the  following: 

From  subroutine  INA77  (Figure  6) , 

(1)  the  input  control  constants  of  namelist  CNTRL 

(2)  the  output  label 

(3)  the  name  card  input  format 

(4)  the  data  input  format 

(5)  the  name  card  information 

From  subroutine  INB77  (Figure  7) , 

(1)  record  number,  subject  number,  and  the  first  NV  data 

(2)  values  for  the  first  ten  subjects  plus  every  K6th 
(NV  being  the  number  of  variables  to  be  processed 
and  K6  as  specified  by  the  user  or  default  of  100) 

(3)  subject  number  of  last  subject  accepted  for  processing 

(4)  number  of  subjects  accepted  for  processing 

(5)  number  of  subjects  read  in 

From  subroutine  TIPAGE  we  get  the  output  which  immediately 
precedes  the  main  XVAL  printout  shown  previously  as  Figure  1, 

• a title  page  (Figure  8) 

• a table  of  contents  (Figure  9) 

Several  elements  of  the  output  reflect  the  effort  made  in  the 
development  of  the  program  to  provide  signals  to  indicate  the 
existence  of  abnormal  outlying  values.  Clues  to  the  "screwball" 
values  can  be  detected  in  several  places. 

The  listing  of  the  values  themselves  will  often  reveal  the 
most  badly  recorded  or  mispunched  data.  Since  it  is  highly 
unlikely,  if  not  impossible,  that  the  youngest  subject  in  this 
survey  of  military  personnel  was  13^  years  old  or  that  the  heaviest 
subject  weighed  in  excess  of  320  pounds,  these  data  suggest  the 
need  for  (1)  checking  back  to  an  original  source  (such  as  data 
blanks),  (2)  further  editing  techniques  or  (3)  elimination  of 
those  values . 

Aberrant  values  can  be  spotted  in  large  differences  between 
the  computed  and  estimated  standard  deviations  which  are  suggestive 
of  abnormal  values  among  the  extreme  20  values.  Gross  differences 
between  BOT  and  TOP  warrant  examination  of  the  smallest  values 
(if  BOT  exceeds  TOP)  or  the  largest  ones  (if  the  reverse  is  true) . 

We  have  found  Veta-II,  the  measure  of  kurtosis,  to  be  the  most 
sensitive  and  most  useful  for  data  having  a more-or-less  normal 
distribution.  For  such  distributions,  the  measure  of  kurtosis  has 
a value  of  about  3.  A single  abnormal  outlyer  can  raise  this  value 
considerably,  on  occasion  almost  to  its  theoretical  maximum,  N, 
the  sample  size.  On  Figure  1,  notice  the  effect  on  kurtosis  of 
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12/21/77  **  • — XVAL  TEST 

Figure  9.  A table  of  contents 


the  deviant  values  for  variable  1 (135.0),  2 (710.0),  4(1389.0), 
7 (-215.0)  and  8 (198.0  and  251.2).* 


PROGRAM  CONSTRAINTS 

The  permissible  size  of  the  input  data  is  not  limited  by 
the  computational  processes,  but  the  use  of  a somewhat  compact 
output  format  does  limit  the  range  of  values  which  may  be  printed 
out  while  preparing  a reasonably  easy-to-read  output.  The  re- 
sults are,  in  the  main,  listed  with  up  to  six  digits  to  the  left 
of  the  decimal  point.  However,  it  will  generally  make  sense  to 
limit  the  data  to  four  digits  to  the  left  of  the  decimal  point  by 
the  use  of  an  appropriate  input  format.  Subject  (record)  numbers 
are  limited  by  the  output  to  integers  with  a maximum  of  five 
digits . 

The  maximum  number  of  variables  is  specified  by  the  user  by 
entering  this  maximum  in  the  "main  program"  as  the  first  dimen- 
sion of  the  arrays  in  the  first  dimension  statement  (not  to 
exceed  the  specified  dimensions  of  arrays  X,  NAY,  and  A,  which 
currently  limit  the  total  number  of  variables  to  205) . The 
number  of  subjects  (or  data  records)  is  not  limited. 

XVAL  does  not  inform  its  user  that  a data  point  is  wrong, 
only  that  it  does  not  fit  the  distribution.  No  alternative 
value  is  suggested  and  therefore  some  knowledge  of  the  data 
and/or  more  intensive  data  analysis  is  often  required  for 
proper  editing. 


SOME  NOTES 

The  following  notes  may  be  of  some  value  in  understanding 
how  the  program  functions. 

Duplicate  variable  names:  The  routine  which  prepares  the 
table  of  contents  will  print  out  any  instances  of  duplicate 
variable  names. 

Calculation  of  summary  statistics:  The  data  for  any  variable 
are  stored  until  there  are  20  of  them.  The  mean  of  these  20 
values  is  computed,  rounded  to  an  integer,  and  designated  as 
A (I, 3).  The  four  summations, 

S ( I , K)  = E(X(I)  - A (I , 3) ) K for  K=l,  2,  3,  and  4 
are  computed.  Finally, 


* Age  in  tenths  of  years,  weight  in  tenths  of  kilograms,  stature 
in  millimeters. 
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A=S(I,1)/N,  B=S(I,2)/N,  C=S(I,3)/N,  D=S(I,4)/N 


and , 


Mean  = A + S (1 , 3) 

Standard  Deviation  = /B  - A*A 

Veta  I = (C  - 3 BA  + 2A3)/(Std  Dev) 3 

Veta  II  = (D  - 4CA  + 6 BA 2 - 3A4)/(Std  Dev) 4 

Estimates  of  mean  and  standard  deviation:  For  any  symmetric 
distribution — normal  or  otherwise — the  mean  of  a set  of  data 
truncated  equally  at  both  ends  is  itself  an  unbiased  estimate 
of  the  mean  of  the  untruncated  set.  Thus,  the  mean  of  the  N-20 
values  is  the  reported  estimate  of  the  mean. 

Any  truncation,  however,  will  result  in  a set  of  data  with 
less  variability  than  the  original  one.  For  a distribution 
of  known  form,  the  effect  of  this  truncation  can  be  estimated. 
Here,  we  assume  a normal  distribution,  and  "correct"  the  actual 
deviation  computed  using  the  N-20  values  by  dividing  it  by 


^1.0124-62.57892/ (N+100) -2 . 57827//N 

Determination  of  interval  widths  and  "minimum"  values: 
Several  programs  which  are  often  used  for  further  processing  of 
the  data  prepare  univariate  or  bivariate  frequency  tables.  These 
programs  require  information  as  to  appropriate  starting  points 
and  interval  widths;  this  information  is  provided  by  the  values 
supplied  by  XVAL  on  the  name-range  cards. 

Several  policy  decisions  are  involved  in  the  way  XVAL  com- 
putes these  values.  First,  the  maximum  of  50  intervals  for  uni- 
variate and  30  intervals  for  bivariate  tables  were  chosen  on 
the  basis  of  being,  in  general,  statistically  adequate  and  about 
as  many  as  were  practical  if  readable  printouts  of  the  ultimate 
tables  were  to  be  obtained.  (Other  maximums  can,  of  course,  be 
specified  on  the  namelist  card.) 

Secondly,  "reasonable"  interval  widths--say,  20.0  instead  of 
18.72 — were  presumed  to  be  desirable.  Consistent  with  the  range 
of  most  of  the  data  for  which  the  program  was  developed,  15 
possible  interval  widths  are  listed  in  XVAL  as  array  WYD  0.1,0. 2, 
0 . 3 , 0 . 5 , 1 . 0 , 2 . 0 , 3 . 0 , 5 . 0 , 10 . 0 , 15 . 0 , 20 . , 25 . , 30 . , 50 . , 80  and  the  pro- 
gram selects  the  smallest  of  these  values  which  will  not  require 
more  than  the  specified  maximum  number  of  intervals.  The  first 
four  (fractional)  values  are  not  considered  appropriate  unless 
the  value  of  IP  in  namelist  XVAL  is  otherwise  specified.  If  the 
largest  of  these  values  is  not  large  enough,  an  interval  width 
which  is  an  integer  times  100  is  used. 
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The  present  selection  of  possible  interval  widths  can  be 
easily  changed  by  replacing  the  array  WYD. 

Third,  to  minimize  the  effects  of  certain  types  of  round^ 
ing  errors  and  to  cause  tables  with  equal  interval  widths  and 
overlapping  ranges  to  coincide  in  their  common  range,  the  follow- 
ing conventions  have  been  used  in  computing  the  lower  limit  of 
each  interval  for  the  univariate  tables,  W being  the  interval 
width  in  computing  a maximum  50  intervals,  and  K an  integer: 

• W less  than  1.0,  K = 0.05 

• W = 1,  2,  or  3,  K = 0.5 

• W = 5 or  more,  K = 2.5 

The  lower  limit  of  the  first  interval  is  chosen  as  the  appropri- 
ate value  of  K,  below  an  integer  multiple  of  the  interval  width, 
W.  Thus,  for  example,  in  Figure  2 we  see  that  for  "cervicale 
height"  with  the  range  1347-1675,  the  interval  width  is: 

1675.0  - 1347.0  „ 


gives  W = 10.0 

134*10  = 1340.0 


is  the  largest  multiple  of  10  that  is  less  than  or  equal  to 
1347.0;  thus  the  chosen  bottom  of  the  first  interval  is: 

134*W  - K = (A ( I , 1 ) ) 

1340.0  - 2.5  = 1337.5 


Occasionally,  this  method  of  establishing  the  intervals  will 
require  a wider  interval  than  would  otherwise  be  necessary.  For 
example,  the  range  135-484  could  be  covered  by  50  intervals 
five  units  wide  (135-139,  140-144,  etc.)  but  XVAL's  conventions 
call  instead  for  intervals  10  units  wide  (132.5-141.5,  142.5- 
151.5,  etc.)  as  follows: 


484.0  - 135.0 
50 


4.98  gives  W = 5 


27*5.  = 135. 


135.  - 2.5  = 132.5 


as  the  chosen  bottom  of  the  first  interval.  However,  the 
range  484.0-132.5  cannot  be  covered  with  50  intervals  of  5., 
thus  the  choice  of  10.  for  the  interval  width. 

Intervals  for  bivariate  tables  are  based  on  the  same  initial 
lower-limits  as  those  determined  for  the  univariate  tables. 
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The  "minimum"  values  as  punched  on  the  name-range  cards  are 
actually  these  lower  limits  and,  as  such,  will  almost  always  be 
less  than  the  actual  minimums  of  the  data.  The  values  punched  as 
the  maximums,  on  the  other  hand,  will  be  true  maximums. 


USE  OF  THE  PROGRAM 

The  "main  program"  does  all  the  statistical  analysis  and 
printing  of  results.  Normally,  the  only  changes  that  would  be 
made  would  be  to  the  dimensions  of  arrays  Z,  NSN,  MISS,  and  S. 
The  first  dimensions  of  these  arrays  should  all  be  greater  than 
or  equal  to  the  maximum  number  of  variables  to  be  processed. 

(As  previously  noted,  program  XVAL,  in  its  present  state,  limits 
the  total  number  of  variables  to  205.) 

The  program  uses  a BLOCK  DATA  subroutine,  a title  page  and 
table  of  contents  subroutine  (TIPAGE) , and  the  INA77  and  INB77 
input  routines,  as  well  as  utility  subroutine  NUNU99  (see  Figure 
10  for  the  deck  setup) . The  contents  of  common  blocks  /DATUMS/ 
and  /HEAD/,  except  for  constants  NMAX  and  NPG  and  array  A,  must 
be  supplied  to  the  main  program  before  the  data  processing  can 
take  place.  In  its  present  state,  INA77  and  INB77  supply  these 
values  as  follows: 

From  INA77  we  get, 

• WHEN  - the  date 

• HDG  - an  80-character  label 

• NV  - the  total  number  of  variables  to  be  processed 

• NAY  - the  array  containing  the  18-character  vari- 

able names 

From  INB77  we  get,  each  time  it  is  called, 

• NSUB  - the  subject  or  record  number 

• X - the  array  containing  the  NV  data  values 

for  processing 

Descriptions  of  these  routines  begin  on  page  40  for  INA77 
and  on  page  49  for  INB77. 

The  Input 

The  data  stream  is  normally  in  the  following  form:  called 
from  the  "main  program" — 

• the  namelist  XVAL 

• the  unit  cards,  if  UNIT  was  set  / 0 in  namelist  XVAL 
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called  from  INA77-- 


• the  namelist  CNTRL 

• a heading  card  for  the  survey 

• the  date,  if  IWHEN  set  ^ 0 in  namelist  CNTRL 

• a format  for  the  name  cards  if  Nl  set  ^ 0 in  CNTRL 

• N2  cards  with  the  format  for  the  data 

• finally,  the  name  cards,  the  last  one  for  variable 

LN  as  set  in  namelist  CNTRL 

called  from  INB77 — 

• for  each  data  set  a subject/record  number  and  NW 
data  values  read  from  unit  NT 

• possibly  a pseudo-data  record  with  a negative  sub- 
ject number,  or  an  end-of-file  card 

The  general  data  deck  layout  is  shown  in  Figure  11. 

The  following  comments  apply  to  these  items: 

(1)  The  namelist  XVAL  and  its  default  values: 

ML  - if  ^ 0 missing  values  are  listed.  If  ML  >0,  up 

to  ML  lines  will  be  printed.  If  ML  <0,  there  is  no 

limit  (default:  0). 

MS  - if  >0,  non-consecutive  subject  numbers  are 
listed  (default:  0). 

NR  - if  >0,  name-range  cards  are  punched  (default:  0) . 

NP  - maximum  number  of  intervals  for  univariate 
tables  (default:  50). 

NQ  - maximum  number  of  intervals  for  bivariate 
tables  (default:  30). 

IP  - the  first  value  of  WYD  to  be  considered 
(default : 5) . 

XKST  - the  usual  factor  for  converting  to  metric 
output  (default:  .1). 

YKST  - the  usual  factor  for  converting  from  metric 
to  English  (default:  .3937008). 

NUNIT  - the  number  of  unit  cards  (default:  0) . 

The  format  of  the  namelist  XVAL  is  as  follows  (see  Figure  12a) : 

column  1 - a blank 

column  2 - a $* 


* This  symbol  is  for  use  of  the  program  on  a CDC  computer;  other 
computers  will  have  different  symbols  for  this  purpose. 
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/ 

• 

/ 

9) 

subject  data 

/ 

/ 

8) 

name (-range)  cards 

/ 

• 

/ 

7) 

N2  data  format  cards 

/ 

6) 

name-range  card  format 

/ 

5) 

a date 

/ 

4) 

a run  label 

/ 

3) 

namelist  CNTRL 

/ 

* 

/ 

2) 

NUNIT  unit  cards 

/ 1)  namelist  XVAL 


Figure  11.  The  XVAL  input  deck. 
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columns  3-6  - the  word  XVAL 
column  7 - a blank 

After  column  7 come  none,  all,  or  any  combination  of  the  control 
variables  in  the  form  X=4.2,N=13,  ...  , the  last  one  followed  by 

a $.  Note  that  although  imbedded  blanks  are  acceptable,  on  some 
machines  they  will  be  considered  as  zeros  when  they  occur  between 
a variable  value  and  the  following  comma.  Thus,  a namelist  string 
of  the  form  X=4.2,N=13  ,1=3$  might  be  interpreted  as  X=4 . 2 N=130 

1 = 3. 


Figure  12b  is  the  namelist  XVAL  input  which  contributed  to 
the  previous  nine  figures  as  follows: 

ML=-1  caused  the  unlimited  missing  value  list  of 
Figure  4a  to  be  printed 

MS=-1  caused  the  missing  subjects  list  of  Figure  4b 
to  be  printed 

NR=1  caused  the  punching  of  the  name-range  cards 
of  Figure  3 

NP=45  and  NQ=25  caused  the  values  for  INTV1  and  INTV2 
on  Figures  2 and  3 to  be  calculated  for  a 
maximum  of  45  and  25  intervals,  respectively. 

IP=4  would  have  allowed  fractional  interval  widths 
beginning  with  WYD(4)  if  any  of  the  variables 
had  the  appropriate  range 

NUNIT=4  sets  a switch  to  allow  the  reading  in  of 
four  UNIT  cards  to  vary  the  metric/English 
conversion  factors  for  the  variables  specified 
on  these  cards 

(2)  If  printing  and/or  punching  of  the  proper  conversion 
factors  is  of  some  concern  to  the  programmer,  "unit"  cards 
may  be  needed.  The  program  automatically  assumes  all  values 
are  in  millimeters  and  desired  units  are  centimeters  and 
inches  and  assigns  0.1  (XKST)  and  0.3937008  (YKST)  as  con- 
version factors.  If  this  isn't  true,  two  things  can  be  done: 
XKST  and  YKST  can  be  reassigned  for  all  variables  in  namelist 
XVAL  or  XKST  and  YKST  can  be  changed  for  specific  variables 
by  specifying  the  appropriate  constants  on  a unit  card.  For 
NUNIT  / 0,  NUNIT  unit  cards  are  read  in, containing  on  each 
card  MISS (L,l) ,MISS (L,2) ,S(L,1) , and  S(L,2).  This  says 
that  all  variables  from  variable  MISS(L,1)  to  variable 
MISS (L, 2)  will  have  metric  (or  primary)  conversion  constant 
S (L , 1 ) and  English  (or  secondary)  conversion  constant  S(L,2). 

These  unit  cards  are  of  the  form  (see  Figure  12c) : 

columns  1-5  MISS(L,1) 

columns  6-10  MISS (L, 2) 
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columns  11-20  S(L,1) 

columns  21-30  S(L,2) 

Figure  12d  shows  the  "unit"  cards  used  to  produce  the 
conversion  factors  (CF1  and  CF2)  in  the  summary  statistics 
of  Figure  2 and  the  name-range  cards  of  Figure  3: 

• card  1 would  convert  variable  1,  age,  in  tenths 
of  years  to  years  for  primary  and  secondary  units 

• card  2 would  convert  variable  2,  weight,  in  pounds 
to  kilograms  for  primary  units  then  back  to  pounds 
for  secondary  units 

• card  3 would  leave  variable  3,  grip  strength,  as 
kilograms  for  primary  units  and  convert  to  pounds 
for  secondary  units 

3 . 

• card  4 would  convert  variable  8,  he ight/ /weight , 
to  centimeters/^kilograms  from  millimeters/^pounds 
for  primary  units,  then  to  inches/^pounds  for 
secondary  units 

Variables  4-6  assume  the  values  XKST(O.l)  and  YKST (0 . 3937008) 
converting  millimeters  to  centimeters  primarily,  and  to  inches 
secondarily . 

(3)  The  namelist  CNTRL  contains  various  control  values  for 
inputting  and  processing  the  data  (see  the  INA77  description) . 

(4)  The  heading  read  inis  printed  out  by  XVAL  on  each  page 
of  output. 

(5)  If  a date  function  is  not  available  on  the  system,  or 
a different  date  is  desired  on  the  printout,  a date  can  be 
read  in  here  and  it  appears  along  with  the  heading  on  each 
printed  page.  If  no  date  is  assigned,  the  field  is  filled 
with  blanks. 

(6)  If  no  format  for  the  name  cards  is  read  in,  the  standard 
format  (see  the  INA77  description)  is  assumed. 

(7)  The  data  format  should  read  NSUB  as  an  integer  and  NW 
floating  point  numbers. 

(8)  The  name  cards  at  this  point  should  contain  a variable 
number  and  a maximum-of-18-characters  variable  name. 

(9)  NSUB  is  the  subject/record  number,  and  (X (I) , 1=1 ,NW) 
are  the  input  data. 

Three  ways  are  available  to  specify  the  number  of  data 
records  which  will  be  processed: 
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• The  number  is  given  as  IRR  on  the  namelist  card. 

• The  end  of  the  data  is  signalled  by  the  reading  of 
a pseudo-data  record  with  a non-positive  value  for 
NSUB.  (Note  that  if  each  data  record  consists 

of  H-cards,  the  pseudo  record  must  also.) 

• The  end  of  the  data  is  signalled  by  the  reading 
of  an  end-of-file  card. 


The  choice  is  up  to  the  user. 


The  Non-array  Variables 

Variables  with  preceding *  *'s  refer,  in  effect,  to  the 
I ' th  variable. 

*AWID  - a computed  univariate  - table  interval  width. 
BLANK  - a utility  variable  equal  to  4H*/v^. 


*BWID  - a computed  bivariate  - table  interval  width. 

*DEL  - the  lower  limit  of  the  table  (A (1,1))  is  chosen 

as  the  [largest  multiple  of  AWID]  - DEL  which 
is  less  than  or  equal  to  the  minimum  value.  If 
AWID  is  greater  than  or  equal  to  5,  DEL  = 2.5  to 
minimize  rounding  errors,  and  to  make  interval 
end  points  for  related  variables  coincide.  If 
AWID  is  less  than  1.,  DEL  = .05.  Otherwise, 

DEL  = .5. 

I - a subscript. 

IP  - the  first  value  of  WYD  to  be  considered. 


J,  JJ,  K,  KE,  KK,  KL,  KLM , L,  LL,  M - subscripts. 


ML 


MM 

MS 

NSUB 


if  not  equal  to  0,  missing  values  are  listed. 

If  ML>0 , up  to  ML  lines  (6  subject  numbers  per 
line)  will  be  printed.  If  ML<0,  there  is  no  limit. 

a subscript. 

if  >0,  non-consecutive  subject  numbers  are  listed, 
record  number  of  the  previous  record. 


MX,N  - a subscript. 

NKQ  - a utility  variable  (=  NK(I)-1  for  NK(I)<20). 


NMAX  - the  maximum  number  of  variables  able  to  be  handled 
in  the  current  run.  It  is  determined  by  the  "main 
program"  as  the  first  dimension  of  array  Z. 
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NP 

maximum  acceptable  number  of  intervals  for 
univariate  tables  used  in  computing  AWID. 

NPG 

the  current  output  page  number  (the  table  of 
contents  being  page  1) . 

NQ 

maximum  acceptable  number  of  intervals  for 
bivariate  tables  used  in  computing  BWID. 

NR 

if  not  equal  to  0,  name-range  cards  are  punched. 

NSUB 

the  current  subject  number  as  read  in  from 
the  data. 

NUNIT  - 

the  number  of  unit  cards  to  be  read  in. 

NV 

number  of  variables  being  processed. 

* RANGE  - 

the  range  of  the  data. 

S2 

a utility  variable  (=S (L, 1) **2 ) . 

*UMB 

(n-20) . 

WHY 

a utility  variable,  used  in  computing  EM. 

XKST 

the  usual  constant  for  converting  this  survey's 
data  to  metric  output. 

*XMID  - 

the  range  between  the  tenth  largest  and  tenth 
smallest  data  values. 

XN 

a utility  variable  (=  NK(I)). 

YKST 

the  usual  constant  for  converting  this  survey's 
data  from  the  converted  metric  output  to 

Englist  output. 

1 

a 

tS3 

* 

a utility  variable  (=  RANGE/NP) . 

*Z9 

a utility  variable  (=  Z (I ,L) -A (I , 3) ) . 

The  Array  Variables 

The  dimensioned  variables  are  those  listed  below.  Those 
dimensions  given  as  205  can  be  changed  to  reflect  the  number  of 
variables  to  be  processed. 

A(205,7)  - the  constants  for  the  name-range  card: 


A (1,1)  = lower  limit  of  1st  interval 

A (I,  2)  = maximum 

A (I,  3)  = approximate  mean 
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A (1 , 4 
A ( 1 , 5 
A (1 , 6 
A (I  / 7 

APE  (14) 

BOT(l) 

CST (205,2) 

EM  (1) , ESD  (1) 

HDG (20) 

MISS  (205, 6) 

MQ(1) 

NAY (205,5) 

NK  (1) 

NSN (205 , 20) 

S (205,5) 

TOP (1) 

V(l) 


) = bivariate  table  interval  width 

) = univariate  table  interval  width 

) = metric  conversion  constant 

) = English  conversion  constant 

- storage  array  for  some  output  headings. 

- equivalenced  with  MISS(1,5),  it  contains 
a measure  of  the  variation  within  the 
ten  smallest  values. 

- metric  and  English  conversion  factors 
(equivalenced  with  A(I,6)). 

- equivalenced  with  MISS(1,1)  and  MISS(1,2), 
they  contain  the  mean  of  the  mid  n-20 
values  and  the  standard  estimate  from 
these  values. 

- a suitable  heading  for  the  survey. 

- MISS(I,J)  is  the  record  number  for  the 
K-th  record  containing  a zero  value  for 
X ( I ) (J=K  mod  6).  MISS  actually  uses 
cells  committed  to  EM,  ESD,  V,  TOP,  BOT, 

MQ,  variables  not  used  until  values  of 
MISS  are  listed. 

- the  number  of  stored  values  of  MISS. 

- (NAY ( I , J) , J=l,5)  is  a maximum  18-character 
name  for  variable  I . 

- equivalenced  with  S(l,5),  it  contains 
the  number  of  processed  values  of  X(I). 

- NSN(I,J)  is  the  record  number  associated 
with  Z (I , J) . 

- (S(I,J),  J=l,4)  contains  the  summations 
for  (X-A(I,3))  to  the  J-th  power  for 
variable  I.  Eventually  the  mean,  stand- 
ard deviation,  and  the  measures  of  skewness 
and  kurtosis  are  stored  here.  S(I,5)  con- 
tains the  sample  size  for  variable  I. 

- equivalenced  with  MISS(1,4),  it  contains 
a measure  of  the  variation  within  the  ten 
largest  values. 

- equivalenced  with  MISS(1,3),  it  contains  the 
coefficient  of  variance  for  each  variable. 
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WHEN  (2) 


the  date. 


WYD(15)  - a list  of  acceptable  interval  widths. 

X(205)  - the  basic  data  for  each  subject. 

Z (205,20)  - the  ten  largest  and  ten  smallest  values 

for  each  variable  in  descending  order. 


THE  SUBROUTINES 

As  noted  earlier,  the  subroutines  supply  the  contents  of 
the  common  blocks  /DATUMS/  and  /HEAD/,  except  for  constants 
NMAX  and  NPG  and  array  A.  Also  included  in  the  subroutines  are 
a title  page  and  table  of  contents  (TIPAGE) , input  routines 
INA77  and  INB77,  and  utility  subroutine  NUNU99.  These  sub- 
routines are  used  for  data  initialization  and  some  preliminary 
analysis  of  the  data. 

Subroutine  BLOCK  DATA 


The  BLOCK  DATA  subroutine  provides  default  values  and 
initial  values  to  subroutine  INA77  and  INB77  as  follows: 


BLOCK  DATA 

CDMMDN/HEAD/H0&< 2U) ,NP&,WHEN(2) 

COrtMON/ IN/ FI (20)  , F2  < 10  0 ) » NL { 1 b)  , I SUE, INK 
OATA  ISUE, INK, HHEN/-1, 0,2*4H  / 

OATA  NL/ 2* 0, 77777, 5, 100, 0,1, 777, U,l,4* 0,1492,0/ 


OATA 

OATA 

OATA 

OATA 

OATA 


OATA  F1/4H  < I4,4H,  2X,  , 4HWA4,  ,4HA2,3, 4HF8 .2 , 4H  , 2F6 , 4H  . 2, 2 , hHF  1 0.  , 2HDATA 
* 7 ) , 1 1* 1 H / DATA 

END  OATA 


10 

20 

30 

40 

50 

60 

70 

80 


ISUE 

INK 


in  INB77  it  counts  the  number  of  subject  data 
sets  read  in.  (Initially  = -1) 


in  INB77  it  counts  the  number  of  subject  data 
sets  passed  to  the  main  program.  (Initially  = 0) 


WHEN  - the  date  is  initialized  to  blanks. 


NL  (20)  - this  array  contains  many  of  the  control  constants 
used  by  the  subroutines  as  well  as  the  main  pro- 
gram. Its  members  are  equivalenced  and  (assigned) 
as  follows: 

NL (1)  - not  used. 

NL ( 2 ) = NW  - total  number  of  variables  (0). 

NL(3)  = NS  - total  number  of  subjects  to  be  read 
in  (77777)  . 

NL ( 4 ) = NT  - the  data  unit  from  which  the  subject 

data  are  to  be  taken  (5) . 
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NL(5)  = K6  - after  printing  out  the  first  10 

subjects'  data,  INB77  prints  out 
every  K6th  subject's  data  (100). 

NL(6)  = LN  - the  variable  number  of  the  last 

name  card  to  be  read  in  (0) . 

NL(7)  = LB  - number  of  the  first  variable  to  be 

checked  for  out-of-range  data  by 
INB77  (1). 

NL(8)  = LT  - number  of  the  last  variable  to  be 

checked  for  out-of-range  data  by 
INB77  (0) . 

NL(9)  = Nl  - number  of  cards  to  be  read  in  con- 

taining the  format  of  the  name 
cards  (0) . 

NL(10)  = N2  - number  of  cards  to  be  read  in  con- 
taining the  format  of  the  data  (1) . 

NL(ll)  = NER  - acceptable  number  of  records  with 

one  or  more  out-of-range  values  (0). 

NL ( 12 ) = IER  - code  for  treatment  of  out-of-range 
values  by  INB77  (0)  . 

NL(13)  = IWHEN  - if  not  equal  to  0,  the  date  (WHEN) 

is  read  in  (0) . 

NL(14)  = IRR  - the  total  number  of  subject  data 

records  to  be  processed  (0) . 

NL(15)  = not  used. 

NL(16)  = KEEP  - a constant  used  for  sorting  out 

unwanted  subject  data  records  (1492). 

NL ( 17 ) = NHDG  - not  used. 

Fl(20)  - this  array  contains  the  default  input  format  for  the 

name  cards  (14 , 2X , 4A4 , A2 , 3F8 . 2 , 2F6 . 2 , 2F10 . 7 ) . 

Subroutine  TIPAGE 


This  routine  contributes  three  basic  outputs: 

• a listing  of  duplicate  variable  names,  if  any  (Figure  13) 

• a title  page  (Figure  14) 

• a table  of  contents  (Figure  15) 

As  soon  as  data  initialization  is  completed,  this  routine  begins 
comparing  variable  names  starting  with  the  first  two  names.  Array 
K is  used  to  store  the  alphabetic  order  of  the  variable  names. 

When  the  order  of  the  first  two  variables  is  complete,  variable 
3 is  compared  with  the  first  two  in  alphabetic  order.  Testing  is 
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12/21/77  **  XVAL  TEST  **  PAGE 

Figure  15.  A table  of  contents. 


only  done  until  the  proper  alphabetic  order  is  found.  Array  K 
is  then  reordered  and  comparisons  begin  for  the  next  variable. 

When  the  alphabetizing  of  the  variable  names  is  complete 
(the  duplicate  names  being  printed  as  they  are  encountered  while 
alphabetizing) , a title  page  and  a table  of  contents  are  printed. 
Control  is  then  returned  to  the  main  program. 

TIPAGE  uses  variable  NV  and  arrays  NAY,  HDG,  and  WHEN  to 
complete  its  task  and  returns  the  page  number  (NPG)  of  the  last 
page  of  the  table  of  contents  to  the  main  program. 

The  Non-array  Variables. 

I,  ILESS1  - subscripts. 

II  &I3  - calculated  page  numbers  for  statistics  are 
stored  here  for  printing. 

12  - calculated  page  number  of  the  first  page  of  the 
summary  statistics. 

15  - calculated  page  number  of  the  first  page  of  the 
statistics  following  the  table  of  contents. 

J - a subscript. 

JL  - the  number  of  variable  numbers  to  be  updated 

in  array  K when  another  variable  name  has  been 
added  alphabetically. 

J96  - the  difference  between  the  number  of  printed 

lines  on  the  last  page  of  a table  of  contents. 
This  is  used  to  space  the  labelling  to  the 
bottom  of  the  last  page. 

KG  - the  number  of  variables  to  be  included  on  the 

left  side  of  the  current  table  of  contents  page. 

KK , L - subscripts. 

LK  - variable  number  of  the  variable  name  currently 
being  printed  on  the  right  side  of  the  page. 

M - a subscript 

N - number  of  variables  to  be  printed  on  the  current 
page . 

NM  - a subscript. 

NMAX  - not  used  here. 

NPG  - the  current  page  number. 
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NPP  - the  number  of  variables  per  page  of  regular 
statistical  output. 

NSUB  - not  used  here. 

NV  - the  number  of  variable  names  to  be  processed. 

The  Array  Variables. 

A(205,7)  - not  used  here. 

HDG(20)  - an  80-character  label. 

K(205)  - the  variable  numbers  are  stored  here  in  the 
order  of  the  alphabetic  arrangement  of  the 
names . 

NAY(205,5)  - the  maximum  18-character  variable  names. 

WHEN (2)  - the  date. 

X(205)  - array  K uses  the  space  allocated  to  array  X 
during  the  execution  of  this  routine. 

Subroutine  INA77 


This  routine  brings  in  the  control  constants  used  in  all 
the  other  routines,  as  well  as  formats,  a heading,  and  variable 
name-range  information. 

The  Input. 

This  input  routine  reads  in: 

• a list  of  program  parameters  (Figure  16a) . 

• a header  label  for  use  in  the  ultimate  output  (Figure  16b) . 

• possibly  a date  for  labelling  the  output  (Figure  16c) . 

• possibly  a format  for  reading  in  the  variable  name-range 
information  (Figure  16d) . 

• the  data  format  (Figure  16e) . 

• the  variable  names  and  perhaps  some  range  and  data  con- 
version information  (Figure  16f) . 

The  program  parameters  are  read  in  using  a NAMELIST.  The  first 
card  providing  these  constants  begins  ^$CNTRL.  All  cards  must  be 
blank  in  column  1,  and  the  list  must  end  with  a dollar  sign  (see 
Figure  17) . Between  the  indicated  beginning  and  end,  any  variable 
on  the  NAMELIST  can  be  specified  in  the  form: 

NV  = 13 

Such  specifications  are  separated  by  commas.  (NOTE:  blanks  immed- 
iately following  a constant  value  are  interpreted  by  some  machines 
as  zeros.  This  NV=13^,  would  give  NV  a value  of  130  on  some 
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Figure  16.  The  INA77 
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machines.)  They  do  not  have  to  occur  in  any  special  order,  and 
any  combination  of  the  listed  parameters — from  one  to  all — may 
be  included. 

The  listed  parameters  are  these: 

NV  = the  number  of  variables  to  be  processed. 

Default  value :NW. 

NW  = the  number  of  variables  in  each  data  record. 

Default  value :NV.  (Ordinarily,  it  is  possible 

for  NV  to  equal,  be  greater  than,  or  be  less  than  NW.) 

NS:  number  of  data  records  to  be  read  in.  If  no  value 

for  NS  is  specified,  the  program  will  continue  until 
either  an  end-of-file  is  encountered  or  a pseudo-data 
record  with  a negative  subject  number  is  read.  Actual 
default  values : 77777 . 

K6 : the  input  lists  the  first  and  every  (K6)th  data  record. 

Default  value: 100. 

NT:  the  designated  unit  number  the  actual  subject  data 

values  are  to  be  read  from.  Default  value: 5 (this 
generally  denotes  card  input) . 

LB , LT , NER , IER : these  four  constants  control  the  checking 
that  the  data  lie  within  the  appropriate  ranges. 

• Values  of  X(I)  for  LB<I<LT  are  checked.  If  LT=0, 
no  checking  is  done.  Default  values : LB=1 ,LT=0 . 

• The  program  aborts  if  the  number  of  records  con- 
taining one  or  more  out-of-range  values  exceeds 
NER.  Default  value : 0 . 

• The  action  taken  by  the  program  on  finding  an  out- 
of-range  value  is  specified  by  IER: 

• IER=0 . The  record  containing  the  faulty  value 

is  passed  over. 

• IER=1 . The  faulty  value  is  replaced  by  zero. 

• IER>1.  The  faulty  value  is  replaced  by  the 

approximate  mean  value  (A(I,3)). 

Nl:  if  not  equal  to  zero,  a format  card  for  the  name-range 
information  is  read.  Otherwise,  the  format  (I4,2X,4A4, 
A2 , 3F8 . 2 , 2F6 . 2 , 2F10 . 7 ) is  implied.  Default  valuerO. 

N2 : the  number  of  format  cards  for  reading  the  data. 

Default  value:l.  Maximum=5. 

LN : if  any  name-range  cards  are  to  be  read  in,  LN  is  the 
variable  number  of  the  physically  last  such  card.  (A 
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convenient  value  is  1 - with  the  card  for  variable  #1 
put  at  the  back  of  the  group  of  cards.)  Default  value: 
MAX (NV,NW) . 

IWHEN : on  some  machines  the  date  stored  in  array  WHEN  can 

be  assigned  by  calling  an  intrinsic  function.  However, 
when  a date  other  than  the  current  one  is  desired  or  no 
such  function  exists  (as  on  our  IBM  machine) , a value 
of  IWHEN  not  equal  to  zero  will  trigger  the  program  to 
read  the  date  in.  Array  WHEN  is  initially  set  to 
Hollerith  blanks.  Default  value :0. 

IRR  = the  total  number  of  subject  data  records  to  be 
processed.  Default  value :NS. 

NHDG:  not  used. 

Figure  16a  shows  the  input  namelist  CNTRL  card  which  con- 
tributed to  the  outputs  of  Figures  1-9  as  follows: 

(1)  The  number  of  variables  to  be  processed,  NV, 
equals  8. 

(2)  Nl=l  causes  the  reading  of  the  name  card  format  of 
Figure  16d. 

(3)  Subject  data  printout  will  consist  of  the  first  10 
subjects'  data  plus  every  10th  (K6th)  as  seen  in 
Figure  7. 

(4)  LN=8  implies  that  variable  name-range  cards  will  be 
read  in  until  variable  8 is  encountered,  thus  the 
input  of  Figure  16f. 

(5)  IWHEN=1  calls  for  the  reading  in  of  the  date  of 
Figure  16c. 

(6)  NW=6  says  that  six  pieces  of  data  will  be  read  in 
for  each  subject. 

The  label  card:  Columns  1-80  contain  the  alphabetic  array 
HDG  which  becomes  available  for  various  labelling  uses.  For  max- 
imum aesthetics,  the  label  should  probably  be  centered  in  this 
field. 


A date:  If  IWHEN  is  set  not  equal  to  zero,  a date  is  read 
in  from  the  first  8 columns  of  the  card. 

The  format:  Format  statements  to  be-  read  in  follow  the 
usual  rules  for  format  statements  included  within  a program 
except  that  the  word  FORMAT  does  not  appear,  and  all  80  columns 
are  available  for  use.  If  N1  is  not  equal  to  0,  there  should  be 
one  card  with  the  name-range  card  format,  followed  by  N2  cards  with 
the  data  format. 
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The  name-range  cards:  At  least  one  card  must  be  read  in 
with  a legitimate  variable  number.  The  maximum  variable  number 
is  205.  Name-range  cards  are  read  in  until  a variable  number 
equal  to  constant  LN  is  read.  Any  part  or  all  of  the  following 
information  can  be  included  here  (see  Figure  18) : 


columns 

1 - 

4 

the  variable  number  (14) 

columns 

7 - 

24 

the  variable  name  (4A4,A2) 

columns 

25  - 

32 

a value  slightly  less  than  or 
the  minimum  (F8.2) 

equal  to 

columns 

33  - 

40 

the  maximum  (F8.2) 

columns 

41  - 

48 

the  approximate  mean  value  (F8 

.2) 

columns 

49  - 

54 

a suitable  interval  width  for 
tables  (F6.2) 

univariate 

columns 

55  - 

60 

a suitable  interval  width  for 
tables  (F6.2) 

bivariate 

columns 

61  - 

70 

a constant  to  convert  the  data  to  the 
usual  metric  units/cm,  kg,  years  (F10.7) 

columns 

71  - 

80 

a constant  to  convert  the  data 
usual  metric  units  to  English 
pounds,  years  (F10.7) 

from  the 
units/inches , 

The  XVAL  program,  while  it  generates  these  values,  does  not  use 
them;  they  are,  however,  commonly  used  in  preparing  univariate  and 
bivariate  frequency  tables  and  in  a number  of  other  programs. 

The  column  numbers  mentioned  reflect  the  default  format  for  this 
input  as  well  as  the  actual  format  of  the  optional  name-range 
card  output  of  the  main  program.  However,  this  input  is  accepted 
by  INA77  in  any  form  as  long  as  a suitable  format  statement  is 
also  read  in  (see  Figures  16d  and  16f) . 

The  Non-array  Variables. 

IER  - initiated  for  use  by  INB77. 

INK,  IRR,  ISUE  - initiated  for  use  by  INB77. 

I WHEN  - if  set  not  equal  to  zero  in  the  namelist,  INA77 
reads  in  a date. 

J - a subscript. 

KEEP,  K6  - initiated  for  use  by  INB77. 

L - a subscript. 

LB  - initiated  for  use  by  INB77. 

LN  - number  of  the  physically  last  name-range  card. 

LT  - initiated  for  use  by  INB77. 

NER  - initiated  for  use  by  INB77. 
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NF2  - equals  20  times  N2  and  represents  the  number  of 
words  of  the  data  format  read  in. 


NHDG  - 

not  used. 

NMAX  - 

the  maximum  allowable  number  of  variables. 

NPG 

not  used. 

NS 

initiated  for  use  by  INB77. 

NSUB  - 

not  used. 

NT 

initiated  for  use  by  INB77. 

NV 

the  number  of  variables  to  be  processed. 

NW 

initiated  for  use  by  INB77. 

Nl 

if  Nl  is  set  not  equal  to  zero,  a format  for  reading 
in  the  name-range  cards  is  read  in. 

N2 

the  number  of  format  cards  to  be  read  in  for  use 
in  reading  the  subject  data. 

The  Array  Variables. 

The  dimensioned  variables  are  those  listed  below.  Those 
dimensions  given  as  205  can  be  changed  to  reflect  the  number  of 
variables  to  be  processed. 

A (205, 7)  - the  constants  which  might  appear  on  the  name- 


range  cards: 

A ( 1 , 1 ) = lower  limit  of  1st  interval, 

A ( I , 2 ) = max imum , 

A (I,  3)  = approximate  mean, 

A (1, 4)  = bivariate  table  interval  width, 

A (1, 5)  = univariate  table  interval  width, 

A(I,6)  = metric  conversion  constant, 

A (1, 7)  = English  conversion  constant. 

FI (20) 

- the  name-range  card  input  format. 

F2  (20) 

- the  subject  data  input  format. 

HDG (20) 

- a suitable  heading  for  the  survey. 
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NAY(205,5)  - (NAY ( I , J) , J=l,5)  is  a maximum  18-character 

name  for  variable  I. 

NL(20)  - this  array  contains  many  of  the  control  constants 

used  by  the  subroutines  as  well  as  the  main  program 

Its  members  are  equivalenced  and  assigned  as  follows 

NL(1)  = not  used, 

NL(2)  = NW  - total  number  of  variables, 

NL(3)  = NS  - total  number  of  subjects  to  be 
read  in, 

NL(4)  = NT  - the  data  unit  from  which  the  subject 

data  are  to  be  taken, 

NL(5)  = K6  - after  printing  out  the  first  ten 

subjects'  data,  INB77  prints  out 
every  K6th  subject's  data, 

NL(6)  = LN  - the  variable  number  of  the  last 

name  card  to  be  read  in, 

NL(7)  = LB  - number  of  the  first  variable  to  be 

checked  for  out-of-range  data 
by  INB77 , 

NL(8)  = LT  - number  of  the  last  variable  to  be 

checked  for  out-of-range  data  by 
INB77 , 

NL(9)  = Nl  - number  of  cards  to  be  read  in  con- 

taining the  format  of  the  name  cards, 

NL(10)  = N2  - number  of  cards  to  be  read  in  con- 
taining the  format  of  the  data, 

NL(ll)  = NER  - acceptable  number  of  records  with 
one  or  more  out-of-range  values, 

NL ( 12 ) = IER  - code  for  treatment  of  out-of-range 
values  by  INB77, 

NL(13)  = IWHEN  - if  not  equal  to  0,  the  date  (WHEN) 
is  read  in, 

NL(14)  = IRR  - the  total  number  of  subject  data 
records  to  be  processed, 

NL(15)  = KEEP  - a constant  used  for  sorting  out 
unwanted  subject  data  records, 

NL ( 16 ) = NHDG  - not  used. 
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WHEN ( 2 ) 


the  date. 


X(205)  - the  basic  data  for  each  subject  are  read  into 

array  X by  subroutine  INB77. 


Subroutine  INB77 


The  Input. 

This  is  our  basic  data  input  routine  and  expects  data  of  the 
form  NSUB , (X ( I ) , (1=1 ,NW)  , where  NSUB  is  a subject/record  number 
and  X(I)  is  the  value  of  the  I ' th  variable  for  subject  number  NSUB. 
Figure  19  shows  the  form  of  the  data  used  in  generating  Figures  1-4 
and  7 as  read  in  by  the  format  of  Figure  16a.  Note  that  the  last 
card  is  blank;  thus,  NSUB  will  be  read  as  zero  and  reading  of  the 
data  will  terminate. 

Data  Manipulation. 

For  each  data  record,  this  routine  reads  in  the  data  record 
and,  usually,  transmits  the  record  to  the  calling  program  after 
it  does  one  or  more  of  the  following: 

• makes  any  necessary  data  changes  through  a call  to 
subroutine  NUNU99 , 

• the  first  ten  records  plus  every  K6th  are  listed 
as  specified, 

• if  LT  was  assigned  in  the  namelist  XVAL,  each  datum 
from  variable  LB  to  variable  LT  is  checked  to  see  that 
it  falls  within  the  range  specified  on  the  name-range 
cards  read  in  from  subroutine  INA77.  When  a value  outside 
that  range  is  detected  and  the  number  of  allowable  records 
with  out-of-range  values  (NER)  has  not  been  exceeded,  one 
of  three  things  can  happen,  depending  on  the  value  of  IER: 

• if  IER=0 , the  record  is  rejected  and  a new  one  used, 

• if  IER=1,  the  out-of-range  value  is  set  equal 
to  zero, 

• if  IER>1,  the  value  is  set  equal  to  the  approxi- 
mate mean  value  for  that  variable  (A(I,3)). 

Note  that  the  third  step  above  is  taken  only  if  the  relevant  pro- 
gram parameter  is  set.  The  first  step  constitutes  a call  to  sub- 
routine NUNU99  which  routinely  does  nothing  except  return  control 
to  INB77.  It  is  in  this  routine  that  we  fabricate  new  variables 
or  eliminate  data  records  by  inserting  sections  of  program  code 
(for  examples,  see  description  of  NUNU99) . 

The  Non-array  variables. 

I - a subscript. 

IER  - code  for  treatment  of  out-of-range  values. 
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Figure  19.  Subroutine  INB77  input. 
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INK 


- counter  that  keeps  track  of  the  number  of 
subject  data  records  processed. 

IRR  - the  total  number  of  subject  data  records  to 
be  processed. 

ISUE  - counter  that  keeps  track  of  the  number  of 
subject  data  records  read  in. 

I WHEN  - used  by  INA77. 

K - a subscript. 

KEEP  - a constant  equal  to  1492  used  in  conjunction 

with  NUNU99  for  eliminating  entire  sets  of  subject 
data . 

K6  - after  printing  out  the  first  ten  subjects'  data, 
INB77  prints  out  every  K6th  subject's  data. 

L - a subscript. 

LB  - number  of  the  first  variable  to  be  checked  for 
out-of-range  data. 

LN  - used  by  INA77. 

LT  - number  of  the  last  variable  to  be  checked  for 
out-of-range  data. 

MSUB  - record  number  of  the  previous  record. 

MX  - utility  variable. 

NER  - acceptable  number  of  out-of-range  values. 

NHDG  - not  used. 

NMAX  - the  maximum  number  of  variables  which  can  be 
handled  in  the  current  run.  It  is  determined 
by  the  "main  program"  as  the  first  dimension 
of  array  Z. 

NPG  - used  by  XVAL  and  TIPAGE. 

NS  - total  number  of  subjects  to  be  read  in. 

NSUB  - subject  number  of  the  current  data  record. 

NT  - the  data  unit  from  which  the  subject  data  are 
to  be  taken. 

NV  - the  number  of  variables  to  be  printed  out  by 
INB77,  and  processed  by  XVAL. 
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NW  - the  number  of  pieces  of  data  to  be  read  in. 


N1,N2  - used  by  INA77. 


The  Array  Variables. 


The  dimensioned  variables  are  those  listed  below.  Those 
dimensions  given  as  205  can  be  changed  to  reflect  the  number  of 
variables  to  be  processed. 


A (205, 7)  - the  constants  which  might  appear  on  the  name- 
range  cards : 

A (1,1)  = lower  limit  of  first  interval, 

A (I, 2)  = maximum, 

A (I, 3)  = approximate  mean, 

A (1, 4)  = bivariate  table  interval  width, 

A (1, 5)  = univariate  table  interval  width, 

A (1, 6)  = metric  conversion  constant, 

A (1, 7)  = English  conversion  constant. 


FI  (20) 

F2 (20) 
HDG (20) 
NAY (205,5) 


the  name-range  card  input  format. 

the  subject  data  input  format. 

a suitable  heading  for  the  survey. 

(NAY(I,J),  J=l,5)  is  a maximum  18-character 
name  for  variable  I. 


NL(20)  - this  array  contains  many  of  the  control  constants 

used  by  the  subroutines  as  well  as  the  main  pro- 
gram. Its  members  are  equivalenced  as  follows: 

NL(1)  = not  used, 

NL(2)  = NW  - total  number  of  variables, 

NL(3)  = NS  - total  number  of  subjects  to  be 
read  in, 

NL(4)  = NT  - the  data  unit  from  which  the 

subject  data  are  to  be  taken, 

NL(5)  = K6  - after  printing  out  the  first  ten 

subjects'  data,  INB77  prints  out 
every  K6th  subject's  data, 
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NL  ( 6 ) 


= LN  - the  variable  number  of  the  last 
name  card  to  be  read  in, 

NL(7)  = LB  - number  of  the  first  variable  to 

be  checked  for  out-of-range  data 
by  INB77, 

NL(8)  = LT  - number  of  the  last  variable  to  be 

checked  for  out-of-range  data 
by  INB77, 

NL(9)  = N1  - number  of  cards  to  be  read  in 

containing  the  format  of  the 
name  cards, 

NL(10)  = N2  - number  of  cards  to  read  in  con- 
taining the  format  of  the  data, 

NL(ll)  = NER  - acceptable  number  of  records  with 
one  or  more  out-of-range  values, 

NL ( 12 ) = IER  - code  for  treatment  of  out-of-range 
values  by  INB77, 

NL (13)  = I WHEN  - if  not  equal  to  0,  the  date 
(WHEN)  is  read  in, 

NL(14)  = IRR  - the  total  number  of  subject  data 
records  to  be  processed, 

NL(15)  = KEEP  - a constant  used  for  sorting  out 
unwanted  subject  data  records, 

NL (16)  = NHDG  - not  used. 

WHEN (2)  - the  date. 

X(205)  - the  basic  data  for  each  subject. 

Subroutine  NUNU99 


This  subroutine's  function  is  to  make  changes  in  the  data. 

It  routinely  does  nothing  except  return  control  to  INB77  (see 
Figure  20).  It  is  in  this  routine  that  we  fabricate  new  variables, 
after  the  regular  variables,  or  eliminate  data  records  by  insert- 
ing sections  of  program  code. 

Figure  21a  shows  the  version  of  NUNU99  which  contributed  to 
the  output  of  Figures  1-4  and  7.  Although  only  six  pieces  of  data 
were  read  in  for  each  subject  by  INB77,  NUNU99  has  created  variables 
7 and  8.  Figure  21b  shows  NUNU99  creating  these  same  two  vari- 
ables, plus: 
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SUBROUTINE  NUNU99  (KEEP) 

C...  THIS  SUBROUTINE'S  FUNCTION  IS  FOR  MAKING  CHANGES  IN  THE  OATA.  IT 
C ROUTINELY  DOES  NOTHING  EXCEPT  RETURN  CONTROL  TO  INB77.  IT  IS  IN 
C THIS  ROUTINE  THAT  HE  FABRICATE  NEW  VARIABLES,  AFTER  ThE  REGULAR 
C VARIABLES,  OR  ELIMINATE  OATA  RECORDS  BY  INSERTING  SECTIONS  OF 
C PROGRAM  CODE. 

C 


C...  WHEN  MAKING  AN  XVAL  RUN,  THIS  ROUTINE  SHOULO  ALWAYS  BE  CHECKED  TO 
C INSURE  NO  UNWANTED  DATA  CHANGES  FROM  A PREVIOUS  RUN  REMAIN. 

n _____ 

w 

c 

CQMMON/QATUMS/X(205) , N A Y ( 20 5 , 5) t A ( 20  5 t 7)  , N V, NMAX , NSUB 


IF  KEEP  .NE.1492,  THE  RECORD  JUST  READ  IS  IGNORED 

KEEP=1492 

RETURN 

ENU 


NUNU  10 
NUN J 20 
NUNU  30 
NUN J 40 
NUNU  50 
NUN J oO 
NUNU  70 
NUNU  60 
NUN  J 90 
NUN J 100 
NUN J 110 
NUN  J 120 
NUN  J 130 
NUN J 140 
NUNU  150 
NUN  J 150 
NUNU  170 
NUN  J 160 
NUN J 190 
NUN  J 200 


Figure  20.  Subroutine  NUNU99 . 
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O O O OOOOOOOOOOOO  O t j O 


SUBROUTINE  NUNU99  (KEEP) 


...  THIS  SUBROUTINE'S  FUNCTION  IS  FOR  MAKING  CHANCES  IN  THE  DATA.  IT 
ROUTINELY  OOES  NOTHING  EXCEPT  RETURN  CONTROL  TO  INB77.  IT  IS  IN 
THIS  ROUTINE  THAT  He.  FABRICATE  NEW  VARIABLES,  AFTER  THE  REGULAR 
VARIABLES,  OR  ELIMINATE  OATA  RECORDS  BY  INSERTING  SECTIONS  OF 
PROGRAM  COOE. 


...  WHEN  MAKING  AN  XV Al  RUN,  THIS  ROUTINE  SHOULO  ALWAYS  BE  CHECKED  TJ 
INSURc.  NO  UNWANTED  OATA  CHANGES  FROM  A PREVIOUS  RJN  REMAIN. 


COMMON/OATUMS/X  (CO 5 > , N A Y ( 205 , 5> , A < 2 05 , 7)  , NV,  NMAX  , NSUB 

IF  <ElP  .NE.1492,  T He  RECORD  JUST  READ  IS  IGNORED 

Ke£P  = l-»92 
X (7)  =X (4) -X  (5) 

IF (X (4) . EQ.  0. .OR. X (5) . ED. 0. ) X (7) =0. 

IF  (X  (4)  .EQ.  0.  .OR. X(2)  . EQ.  0.)X(d)=0. 

IF (X (4) .EQ. 0. .OR.X (2) . EQ. 0 .) RETURN 
X ( 8)  =X(4)/X(2)**(l./3) 

RETURN 

END 


NUN J .0 
NUN J 20 
NUNU  30 
NUN  J 40 
NUNU  50 
NUN J bO 
NUN J 70 
NUNU  80 
NUN J 90 
NUN J 100 
NUN J 110 
NUNU  120 
NUN J 130 
NUN J 140 
NUN J ISO 
NUNU  160 
NUNU  1 7 U 
NUN J 180 


NUN  J 190 
NUNU  2 J 0 


a . 


SUBROUTINE  NUNU99  (KEEP) 


...  THIS  SUBROUTINE'S  FUNCTION  IS  FOR  MAKING  CHANGES  IN  THE  DATA.  IT 
ROUTINELY  DOES  NOTHING  EXCEPT  RETURN  CONTROL  TO  INB77.  IT  IS  IN 
THIS  ROUTINE  THAT  WE  FABRICATE  NEW  VARIABLES,  AFTER  THE  REGULAR 
VARIABLES,  OR  ELIMINATE  OATA  RECORDS  BY  INSERTING  SECTIONS  OF 
PROGRAM  COOE. 


...  WHEN  MAKING  AN  XVA„  RUN,  THIS  ROUTINE  SHOULO  ALWAYS  BE  CHECKED  TO 
INSURE  NO  UNWANTED  DATA  CHANGES  FROM  A PREVIOUS  RUN  REMAIN. 


COMMON/OATUMS/X (205) ,NAY(205,5),A(205,7) , N V , N M A X , NSU B 
IF  KEEP  .NE.1492,  THE  RECORO  JUST  READ  IS  IGNORED 
KEEP=1492 

IF(NSUB.EQ.145)KEEP=1 
IF(NSUB.EQ.16)X(2) =0 
IF ( NSUB • EQ. 12  0 ) X (1) =0 
IF ( NSUB . EQ. 82) X(4>  =0. 

X ( 7 ) =X(4)-X(5) 

IF (X (4) .EQ.O . .OR.X (5) . EQ. 0.) X (7) =0. 

IF (X (4) .EQ.G. . OR.X (2) .EQ. 0. ) X (8) =0. 

IF (X (4) .EQ.O. .OR.X (2) .EQ.  0. ) RETURN 
X (8) =X (4) /X (2) **  (1. /3) 

RETURN 

END 


NUNU  10 
NUNU  20 
NUNU  30 
NUNU  40 
NUNU  50 
NUNU  60 
NUNU  70 
NUNU  80 
NUNU  90 
NUNU  100 
NUNU  110 
NUNU  120 
NUNU  130 
NUNU  140 
NUNU  150 
NUNU  160 
NUNU  170 
NUNU  180 


NUNU  190 
NUNU  200 


b. 


Figure  21.  Editing  with  NUNU 9 9 . 
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• deleting  subject  number  145  from  the  sample, 

• setting  weight  for  subject  16  equal  to  zero, 

• setting  age  for  subject  120  equal  to  zero, 

• setting  stature  for  subject  82  equal  to  zero. 

In  this  manner,  we  have  edited  out  subject  145,  who  appears  to  be 
substantially  smaller  than  anyone  else,  and  also  edited  out  three 
individual  values  that  seem  to  be  undoubtedly  in  error.  The 
result  is  shown  in  Figure  22. 
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AGE  WEIGHT  GRIP  ST*  STATURE  CERVICAL  ACROMION  STATURE-  HEIGHT/W 

ENGTH  fc  HEIGHT  HEIGHT  CERVICALE  EIGHT**l/3 

VALUE  SBJCT  VALUE  SBJCT  VALUE  SBJCT  VALUE  SBJCT  VALUE  SrfJCT  VALUE  SBJCT  VALUE  SBJCT  VALUE  SBJCT 
1ST  SMALLEST  305.0  95  140.0  128  41.0  154  1643.0  173  1392.0  173  1329.0  60  221.0  77  289.3  106 
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Figure  22.  An  example  of  edited  statistics. 


THE  EDITING  PROGRAM* 


The  EDIT  program  is  designed  to  single  out  from  a mass  of 
normally  distributed  data  those  values  which  appear  to  be  ab- 
normal or  aberrant  when  judged  in  terms  of  other  data  for  the 
same  subject.  The  program's  procedure  is  based  on  the  computa- 
tion of  multiple  regression  equations  for  estimating  a given 
variable  in  terms  of  a pair  of  related  variables,  followed  by  a 
comparison  of  the  "actual"  values  with  their  regression  estimates. 
When  the  differences  between  these  values  exceed  a preassigned 
criterion  (expressed  in  terms  of  the  regression  equation's 
standard  error  of  estimate) , the  program  reports  this  fact, 
together  with  a variety  of  relevant  information.  A decision  can 
then  be  made  as  to  how  the  aberrant  value  should  be  treated.  The 
EDIT  program  consists  of  the  main  computer  program,  a series  of 
subroutines  and  a series  of  control  cards  which  relate  to  the 
input  data  and  how  it  will  be  used  in  the  program.  An  EDIT  com- 
puter deck  setup  is  illustrated  in  Figure  23. 


THE  PROGRAM  INPUT 

In  using  the  program,  one  supplies  the  computer  with  five 
sets  of  information:  the  combinations  of  variables;  the  lists 
of  variables  to  be  used  for  comparison;  the  criteria  of  aberrance; 
the  number  of  iterations;  and  the  number  of  subjects  for  equation 
computation . 

The  Combinations  of  Variables 


Variables,  in  sets  of  three,  are  usually  selected  so  that 
the  multiple  correlation  for  each  one  of  the  three  in  terms  of  the 
other  two  is  fairly  high.  Examples  of  useful  combinations  include: 

• StatuAe , i>houldeA  height,  chei>t  height.  These  variables 
have  the  double  virtue  of  being  of  the  same  anatomical 
type  (heights)  and  of  being  overlapping  measurements. 

• depth,  wali>t  bAeadth,  wali>t  clAcum6eAence.  The 
combination  of  breadth,  depth,  and  circumference  is  tied 
fairly  closely  together  by  geometric  restraints.  These 
three  variables  may  not  be  highly  correlated  individually, 
but  the  multiple  correlations  should  be  high. 

• Steep*  clAcum6eAence,  elbow  cd.Acum6e.Ae.nce.,  ^0-ieaA.m  clAcum- 
6e.Aence.-~ three  measurements  of  the  same  type  made  on  the 
same  part  of  the  body . 

• Hand  length,  { \oot  length,  AtatuAe- - three  "long-bone" 
measurements . 


* A complete  printout  of  the  EDIT  program  is  included  as  Appendix  B. 


58 


Figure  23.  EDIT  deck  setup. 
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For  each  combination,  the  program  computes  three  equations 
and  the  related  standard  errors  of  estimate  for  estimating  each 
variable  in  terms  of  the  other  two.  The  order  in  which  the 
three  variables  are  specified  has  no  significance.  Up  to  125 
variables,  in  a maximum  of  100  combinations,  may  be  used  in  the 
editing  program.  A given  variable  may  be  used  in  any  number 
of  combinations. 

It  is  often  desirable  to  create  variables  to  use  in  the 
editing  program  from  read-in  variables . We  might,  for  example, 
create  the  variable,  stature/^weight , from  the  read-in  variables 
of  stature  and  weight.  To  do  this,  it  is  only  necessary  to  write 
the  appropriate  FORTRAN  statement,  insert  this  in  subroutine 
NUNU99  and  change  the  value  of  the  statement  for  NVC  on  the 
$ED  NAMELIST  card  to  the  number  of  computed  variables.  These 
computed  variables  will  appear  on  the  printout  in  continuing 
numerical  sequence  following  the  read-in  variables. 

If  it  is  desirable  to  compare  a variable  with  only  one 
rather  than  two  other  variables,  this  can  be  done  by  specifying 
a constant  as  the  third  variable  in  the  combination  (see  "Com- 
bination Cards"  for  EDIT  program,  page  64)  . A constant  can  be 
created  by  putting  a FORTRAN  statement  in  subroutine  NUNU99. 

The  Lists 


When  a value  is  flagged  as  being  apparently  aberrant,  it  is 
often  desirable  to  compare  this  value  with  a variety  of  other 
data.  Thus,  if  the  stature-shoulder  height-chest  height  combina- 
tion indicates  that  one  of  these  heights  is  out  of  line  with  the 
other  two,  it  might  be  important  to  examine  all  other  height  data 
for  the  subject  in  question.  If  the  waist  breadth-depth-circum- 
ference combination  signalled  a possibly  discrepant  value,  one 
might  wish  to  check  on  the  subject's  weight  and  his  other  torso 
breadths,  depths,  and  circumferences. 

The  purpose  of  the  lists  is  to  specify  the  additional  vari- 
ables with  which  out-of-line  values  from  a given  combination  can 
be  compared.  Numbered  lists  containing  as  many  as  18  variables 
are  read  into  the  computer;  each  combination  should  have  associ- 
ated with  it  the  number  of  one  of  these  lists.  There  may  be  as 
many  as  20  different  lists.  A variable  can,  of  course,  appear  in 
more  than  one  list,  and  a list  can  be  associated  with  any  number 
of  combinations. 

The  Criterion  of  Aberrance 


A datum  is  treated  as  possibly  aberrant  if  the  difference 
between  the  regression  equation  estimate  and  the  recorded  value 
exceeds  a specified  number  of  standard  errors  of  estimate.  The 
choice  of  this  number  (CK)  is  arbitrary  and  depends  somewhat  on 
the  sample  size,  the  number  of  variables  and  combinations,  the 
state  of  the  data  and  so  forth.  If  no  criterion  value  is  speci- 
fied, the  computer  supplies  a default  value  of  3.5.  The  person 
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using  the  program  for  the  first  time  may  do  well  to  use  this  value 
initially,  subsequently  raising  or  lowering  it  as  seems  desirable. 

With  very  small  samples,  the  selection  of  the  criterion  must 
take  into  account  the  fact  that  no  difference  between  "actual" 
and  regression  values  can  ever  exceed  /N-l  standard  deviations 
(N=number  of  subjects).  With  the  criterion  set  at  4.0,  there 
can  be  no  "aberrant"  values  if  N<_17 . 

Number  of  Iterations 


The  program  reads  in  a specified  number  of  data  records, 
accumulating  the  values,  squares,  and  cross  products  necessary 
to  determine  the  regression  equations,  and  storing  the  data 
as  it  proceeds.  Once  the  regression  equations  have  been  computed, 
the  data  are  reentered  into  the  system,  and  the  comparisons  of 
estimated  values  with  actual  values  are  made.  Whenever  an  out-of- 
line  value  is  sensed  for  an  individual,  the  data  for  the  comparison 
which  flagged  the  out-of-line  value  are  removed  from  the  summations. 
After  all  the  comparisons  have  been  made,  the  equations  are 
recomputed,  and  the  process  of  making  the  comparisons  is  repeated. 

The  reason  for  recomputing  the  equations  and  remaking  the 
comparisons  is  that  if  the  data  contain  even  a few  highly 
erroneous  values,  the  regression  equations  based  on  these  values 
will  be  erroneous  and  the  standard  errors  of  estimate  will,  as 
a rule,  be  quite  high.  Ordinarily,  highly  erroneous  values  will 
be  caught  on  the  first  round  of  comparisons.  By  removing  them  from 
the  summations,  one  gets  more  realistic  and  reliable  equations 
and  smaller,  but  still  significant,  discrepancies  can  be  caught. 

The  computer  will  set  the  number  of  iterations  at  two  unless 
otherwise  instructed.  Usually  this  will  be  satisfactory.  If 
the  data  are  quite  "dirty,"  three  iterations  may  be  desirable. 

When  the  data  on  which  the  equations  are  being  computed  have 
already  been  edited,  a single  iteration  may  suffice. 

Number  of  Subjects  Used  in  Computing  the  Regression  Equations 

Unless  contrary  instructions  are  provided,  the  program  will 
base  the  regression  equations  on  the  data  for  all  the  subjects. 

With  large  numbers  of  subjects,  it  may  be  practical  and  conserving 
of  computer  time  to  base  the  equations  on  some  fraction  of  the 
sample,  say  the  first  500  subjects.  Similarly,  when  editing  a 
set  of  data  consisting  of  already  edited  records  followed  by  a 
number  of  unedited  records,  the  most  satisfactory  procedure  is 
often  to  base  the  equation  on  the  clean  or  previously  edited 
data  and  set  the  number  of  iterations  to  unity. 

To  change  the  program  to  use  only  a limited  set  of  subjects, 
set  the  value  of  NEQ  on  the  $ED  NAMELIST  to  the  appropriate 
number  of  subjects.  When  this  option  is  exercised,  the  program 
goes  through  the  specified  number  of  iterations  for  these  first 
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NEQ= ( ) subjects,  and  then  proceeds  to  read  in  and  check  the  data 
for  all  the  remaining  subjects. 

Number  of  Subjects  in  Data  Set 

As  the  program  is  presently  written,  the  maximum  number  of 
subjects  cannot  exceed  9998.  On  the  rare  occasions  when  the 
subjects  exceed  this  number,  the  program  can  be  modified  by 
changing  the  value  for  NTOTAL  (total  number  of  subjects,  plus 
one)  on  the  $ED  NAMELIST  to  the  appropriate  value.  If,  for 
example,  we  wanted  to  modify  the  program  to  accept  a sample  size 
of  10,001,  the  value  of  NTOTAL  would  be  changed  to  at  least 
10,002,  a value  larger  than  the  size  of  the  sample. 

Size  of  Input  Data 

The  number  of  digits  for  any  datum  value  will  have  little 
effect  on  the  calculations  but  at  several  points  the  printouts 
are  fairly  compact.  Therefore,  data  input  values  should  be 
expressed  with  no  more  than  four  places  to  the  left  and  one  to 
the  right  of  the  decimal  point.  (The  program  was  written  so 
that  a 120-character  printer  would  print  its  output.  When  the 
program  is  used  in  conjunction  with  a 132-character  printer,  the 
F6.1's  in  several  of  the  output  format  statements  could  easily 
be  changed  to  F7.1's  or  F7.2's.)  When  the  values  are  either 
too  big  or  too  small,  they  can  be  adjusted  either  by  a suitable 
format  statement  or  by  addition  of  the  NUNU99  subroutine  to  the 
basic  input  program. 

Major  Programming  Changes 

The  total  number  of  variables  that  can  be  used  is  presently 
set  at  125,  but  it  could  be  much  larger  on  many  of  the  larger 
late-model  computers.  To  increase  this  limit,  it  is  necessary 
to  replace  125  with  the  desired  variable  number  value  in  the 
dimensional  statements  for  X,  NAME  and  A in  COMMON/DATUM  in  the 
main  program  and  in  each  subroutine.  It  should  be  pointed  out, 
however,  that  increasing  the  number  of  variables  uses  computer 
memory  and  the  upper  limit  for  such  variables  is,  therefore, 
finite . 

The  number  of  combinations  can  be  increased  beyond  its  present 
100  by  replacing  this  value  with  the  appropriate  number  in  the 
COMMON/COMP  statements  for  JOB,  S,  SS  and  CSQ.  The  value  of 
MAXCOM  must  also  be  changed  to  match  the  new  combinations  value. 

To  increase  the  number  of  lists,  change  the  dimension  in 
COMMON /COMP  for  LIST  from  the  present  20  to  the  appropriate 
value  and  the  IF  statement  in  line  62  of  the  program  to  the 
same  value. 


62 


THE  PREPARATION  OF  PROGRAM  INPUT  CARDS 


The  namelist,  the  first  card  of  data,  resets  the  values  of 
a number  of  control  variables.  The  initial  values  for  all  the 
namelist  variables  are  written  into  the  program.  If  no  changes 
in  these  values  are  needed,  a blank  namelist  is  used  ( $ED  $) . 

If,  however,  changes  are  required  to  handle  a specific  set  of  data, 
the  following  options  can  be  executed  by  specifying  the  approp- 
riate values  in  the  NAMELIST  card: 

CK  - the  test  value  in  number  of  units  of  standard 
errors  of  estimate  away  from  the  regression 
value  that  will  flag  a datum  for  inspection. 

The  value  is  set  at  3.5  units  unless  changed. 

NREP  - the  number  of  iterations  the  program  goes 

through.  (Set  at  two  in  the  program  unless 
changed . ) 

NEQ  - the  number  of  subjects  used  on  which  to  base 

the  regression  equations.  NEQ  is  set  equal  to 
the  total  sample  unless  a value  is  specified 
in  the  namelist. 

NTOTAL  - a number  larger  than  the  total  number  of 

subjects.  (Initially  set  to  9999  in  the  program.) 

NPUNCH  - if  NPUNCH  is  not  equal  to  zero,  the  variable 
number,  variable  name,  the  mean,  and  the 
standard  deviation  for  each  variable  are 
punched  out.  (Initially  set  to  zero  in 
the  program.) 

NPRINT  - if  NPRINT  is  one  or  more,  infomation  on  iterations 
other  than  the  last  is  printed  out  and  a listing 
of  the  combinations  together  with  the  single  and 
multiple  correlations  are  printed  out.  (Initially 
set  to  zero  in  the  program.) 

NCV  - the  number  of  computed  variables.  (Initially  set 
to  zero  in  the  program.) 

Let  us  assume  that  it  is  necessary  to  edit  a new  data  set  for 
10,000  subjects.  As  this  is  the  first  run  on  these  data,  we 
anticipate  that  there  will  be  a number  of  aberrant  values  to  scan. 
We  want  to  set  the  test  value  at  5 rather  than  3.5  and  the  number 
of  repetitions  at  three.  We  want  the  mean  and  standard  deviation 
for  each  variable  to  be  punched  and  each  iteration  to  be  listed. 

For  editing  purposes,  we  wish  to  base  the  equations  on  only  the 
first  500  records  and  to  create  12  new  variables  from  the  data 
being  read  in.  The  NAMELIST  card  is,  therefore,  modified  as 
shown  in  Figure  24a.  Only  a single  modification  was  made  on  the 
NAMELIST  card  used  in  the  sample  program  which  will  be  described 
later  (see  Figure  24b) . 
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$EB  CK=5»  HREP=3»  MEGN500.  NT0TAL=1DQG1 j MPUNCH=1j  HPRINT=i»  HCV=1£$ 


000000000  0000000  000000  00  0 0000  0 000  00000  000000  00  000  000000000000000 

1 2 3 4 9 6 ; 6 9 K.  1112  13  14  IS  16  If  16  19  20  21  22  23  24  29  26  27  26  29  30  31  32  33  34  39  36  37  36  39  40  41  42  43  44  49  46  47  48  49  SO  SI  S2  S3  94  SS  S6  S7  S6  99  60  61  62  63  64  69  66  67  66  69  70  71  7273  74  75  76  77  76  79  60 

11111111111111111111111111111111  11  111  111111111  111111111  111111  1111111111111 


nnnnnn  nnnnoooooooooooooooooooooooooooooooooooooooooooooooooooooooo  ooooooooooon 

Figure  24a.  Modified  NAMELIST  input  card. 


$EB  NPRINT  = 9 $ 


00000000000  00000000000000000000000000000000000000000000000000000000000000000000 

I 2 3 4 S 6 ; 6 9 1011  12  13  14  IS  16  If  16  19  20  21  22  23  24  29  26  27  26  29  30  31  32  33  34  39  36  37  39  39  40  41  42  43  44  49  46  47  46  49  SO  SI  S2  S3  94  SS  96  S3  96  99  60  61  62  63  64  69  6667  69  69  73  71  72  73  24  7S  76  77  76  79  60 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 111  1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 1 1 1 1 I 


Figure  24b.  Slightly  modified  NAMELIST  input  card. 


Combination  Cards 


The  combination  card,  illustrated  in  Figure  25,  is  divided 
into  six  fields:  columns  5-16,  17-28,  29-40,  41-52,  53-64, 

65-76.  Each  field  contains  the  information  for  one  combination 
in  the  following  pattern: 

5th-7th  columns:  number  of  1st  variable 

8th-10th  columns:  number  of  2nd  variable 
llth-13th  columns:  number  of  3rd  variable 
14th-16th  columns:  number  of  associated  list 

The  first  four  columns  (COMB)  are  used  to  identify  the  card  as 
a combination  card. 

COMB  2 3 4 1 3 4 5 1 6 7 8 1 9 10  11  1 £ 12  13  1 1 7 3 1 


000000000000000000000000000000000000000000000  0000000000000000000000000000000000 

1 2 3 4 S 6 ; 9 9 10  31  32  33  14  IS  16  17  16  19  20  21  22  23  24  2S  26  27  26  29  30  31  32  33  34  3S  36  37  36  39  40  41  42  43  44  4S  46  47  46  49  SO  SI  92  S3  94  SS  S6  97  90  99  60  61  62  63  64  6S  66  67  66  69  70  71  72  73  74  7S  76  77  78  79  60 

111111111111111  11111111111  11111111111  1111  11  11  1111  II  111  II  11111111  1111 

Figure  25.  COMBination  card. 

There  is  no  particular  significance  to  the  order  in  which 
the  three  variables  that  make  up  a combination  are  listed,  nor, 
with  one  exception,  is  there  any  significance  to  the  order  of 
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the  combinations.  An  individual  card  can  contain  anywhere  from 
one  to  six  combinations,  using  any  of  the  six  fields  on  the  card. 

The  one  instance  where  the  order  in  which  the  combinations 
are  read  into  the  computer  makes  a difference  is  related  to  the 
table  for  the  standard-score  values  which  precedes  the  list  of 
errors.  This  table  is  based,  for  each  variable,  on  the  mean  and 
standard  deviation  as  calculated  in  connection  with  the  last 
combination  in  which  that  variable  appeared.  If  a variable 
is  missing  from  a substantial  number  of  data  records,  it  may 
be  well  to  put  combinations  involving  this  variable  near  the 
beginning  as  the  standard-score  values  for  variables  grouped 
with  a missing  variable  will  also  be  based  on  the  reduced  sample. 

The  maximum  number  of  combinations  is  set  at  100  by  the 
program. 

List  Cards 

The  list  cards  take  the  form  shown  in  Figure  26: 

columns  1-4  - LIST 

columns  5-7  - number  of  the  list 

columns  8-10  - number  of  first  variable  on  list 

columns  11-13  - number  of  second  variable  on  list,  and  so  on. 

l i ST  1 1 2 3 45678  9 10  11  ic  U 


0 U 00000000000000000000000000000000  0000000000000000000000000000000000000000000 

1 2 J « 5 6 1 3 9 10  II  12  II  II  lb  IS  II  II  >3  20  21  22  21  2«  35  26  21  23  29  30  21  02  00  ol  15  16  01  08  09  Mil  12  1 0 U 15  16  I!  18  19  50  51  52  50  51  55  56  51  58  59  80  61  62  51  61  65  66  61  68  69  19  01  1?  10  01  15  !»  11  M 19  80 

I 11  11  I 11  1111111111111111111111111  11  1 II  1 1 1 1 1 I 1 1 1 1 1 1 1 1 1 | | 1 | 1 1 1 | 1 1 1 1 1 1 1 | 1 M 1 

Figure  26.  LIST  card. 

There  may  be  anywhere  from  1 to  18  variables  on  a list; 
there  may  be  blanks  in  the  arrangement  on  the  cards.  List 
numbers  need  not  be  consecutive  although  the  highest  permissible 
number  is  20. 

The  list  and  combination  cards  are  placed  after  the  $ED 
$ NAMELIST  and  are  followed  by  a blank  card  (see  Figure  27)  . The 
program  checks  the  first  four  columns  for  COMB,  LIST  or  blank;  if 
anything  else  is  found,  the  program  prints  "something  is  wrong 
with  this  combo-list  card."  After  the  program  finds  a blank 
card,  the  next  card  must  be  a $CNTRL  NAMELIST. 

Data  Input 

The  data  input  for  the  EDITing  program  is  identical  to  that 
previously  described  for  the  XVAL  program  (see  sections  on 
INA77,  INB77  and  NUNU99).  A typical  card  setup  begins  with  a 
namelist  CNTRL  card  on  which  we  enter  the  number  of  variables. 
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The  next  card  contains  a heading  to  identify  the  subject  group, 
in  this  case  "USAF  1967  first  150  subjects."  The  next  card  is 
the  format  for  the  data  to  be  punched  on  the  card,  in  this  case 
"(14,  13F5.0)."  Following  this  is  a group  of  cards  each  con- 
taining the  name  and  number  of  individual  variables.  This  group 
is  followed  by  cards  listing  the  individual  subject  data  in 
the  specified  format. 

The  input  setup  described  above  is  shown  in  Figure  27.  This 
setup  is  for  the  CDC  6600  and  will  be  somewhat  different  for 
other  computers. 


USING  THE  EDIT  PRINTOUT 

To  illustrate  the  use  of  the  program  output  for  editing  anthro- 
pometric data,  a sample  printout  has  been  prepared  and  is  shown  in 
Figures  28-31. 

The  first  component  (Figure  28)  shows:  (1)  the  NAMELIST  $ED 

values;  (2)  the  program  parameters  including  the  combinations 
and  the  lists;  (3)  values  for  NAMELIST  $CNTRL;  (4)  identification 
of  the  subject  group  and  listing  of  the  format  for  the  name  cards 
and  the  data  cards;  (5)  list  of  name  cards;  and  (6)  data  from  the 
first  10  and  every  100th  subject. 

Following  verification  of  the  input  as  described  above  and 
shown  in  Figure  28,  the  editing  begins.  The  summary  statistics  for 
the  variables  by  editing  combinations  for  the  first  data  run  are 
listed  (Figure  29).  These  include,  for  each  variable,  the  mean, 
standard  deviation,  the  multiple  correlation  of  the  variable  with 
the  other  two  variables  in  the  combinations  as  well  as  the  regres- 
sion equations,  the  simple  correlations  and  the  sample  size  (N) . 

Figure  30  shows  a listing  of  subjects  that  were  flagged  as 
having  aberrant  values  based  on  the  value  of  CK.  Each 
entry  in  the  list  consists  of  the  subject  number  followed  by 
the  variable  number  and  name,  the  recorded  value,  the  regression 
value,  and  the  number  of  units  of  standard  error  of  estimate  the 
values  differ  for  each  of  the  three  variables  in  the  combination. 

The  summary  statistics  are  again  given  as  computed  after  the 
aberrant  values  flagged  in  the  initial  run  have  been  deleted 
(Figure  31) . Note  that  the  multiple  correlation  coefficients 
have  been  improved  with  the  removal  of  the  aberrant  values. 

With  the  printout  in  hand  the  NAMELIST  values  given  at  the 
head  of  the  printout  are  checked  to  assure  that  the  intended 
input  parameters  were  indeed  used.  The  combinations  and  lists 
are  also  scanned  in  order  to  determine  that  each  of  the  variables 
of  interest  are  included  in  at  least  one  combination  and  in  the 
list  called  for  by  that  particular  combination  (Figure  28,  (2)). 

The  summary  statistics  and  list  of  aberrant  values  can  initially 
be  ignored.  The  second  set  of  summary  statistics  is  scanned 
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Figure  27.  Typical  card  setup. 
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Figure  28.  Input  verification. 


US  AF  1967  FIRST  150  SUBJECTS 
<I4,2X,4A4, A2,3F8. 2, 2F 6.2, 2Fl 0.7) 
(14, 13F5.0) 


n 

<x 

UJ 

or 

o 

oc 

o 

o 

UJ 

oc 

X 


000'300000000< 

oooooooooooo< 
no'J'^oooooooo* 
ooeoooooonoo  ( 
oooooooooooo< 


I i I I I i I I I I l I ( 


* 

o 

o 


UN 


OOOOOO^OOOOoO 
OOOOOOOO  “30000 

ooooooooooooo  IU 

000000-30  30000  X 

O 30000000300  3 ►— 

OO  30000000000  • • • • . • • • • • • tO 

III  I I I I I I • I I I O CP  00  CVJ  « K CO  to  N.  O N.  4 

x0  4DP-.«)CPv0P.<0c0®X 
PJPJPJPJPJPJPJiVJPJPJPj 

OOOoOOOOOOOOO  to 

0000  300  00  3 000  *-« 

• • • • * • • • • • .X 

00300003000  0 0 <0  PO  P.  xO  O O <X>  H qo  N.  rp  t— 

1 I I I I I I I I I I I I <0  CP  « <0  O <H  «o  CP  O CP  CP 

CO 

N- 

• ••  ••  • ••••  * v-4 

tO(VlM«N>«O)OlA0OH 

OOOOOOOOOOOOO  JO  O P.  P>»  O O N-  CO  p*.  CP  p*. 

3000000000000  cvj  <\j  (\J  pj  ro  ro  pj  pj  cvj  pj  rj  • 

o 

OOOoOOOOOOOOO  T. 

I I I I I I I I I I I I I • * • » • • • • • • • 

—4  JV  CP  PJ  <0  OO  ro  ro  IP  O H- 
,D*.vOp.CPOP*«OCOcOxOO 

CVJ  <\J  PJ  CVJ  PJ  fO  CVJ  PJ  PJ  CVJ  PJ  UJ 

"3 

CD 

OOOOOOOOOOOOO  • • • * • • • • • • a "3 

& vO  P»  vD  O 4 xH  4-  CVJ  CO  O (/) 

ooooooooooooo  co  cvj  CP  cvj  ro  4 o cvj  «-i  p- 

I I I I I 1 1 1 I I I I I cvj  ro  cvj  ro  «o  ro  ro  ro  ro  ro  cvj  QC 

o 

Ll. 

U>x-4UVrOOLCVIS.UV(PvO®0 

CP4ouvcvj^r».cvjcPx4irva 

<pe040CPOe0p»a3<0CP«0O 

OOOOOOOOOOOOO  h O 

UJ 

00000C30000000  • • • • • • • • • • • oc 

1 1 1 I 1 I 1 1 1 1 1 1 1 s <h  ro  h s CP  P»  cp  x>  CVJ  O 

roo<P'X>u\op.copj  4 cvj  u. 
rorocvirororocvjcvjrororoo 

O 



CPUV4N-vOt-IcOIJVCP-I4»-i 

OOOOOOOOOOOOO  <D  CVJ  O 4 LTV  xH  CP  CVJ  xH  UV  T^O 

CVJ  CVJ  CVJ  CVJ  CVJ  CVJ  X-4  PJ  pj  PJ  CVJ<Z 

OOOOOOOOOOOOO  UJ 

I I I ((  I I I I I I I I DC 

4 'JV  Pj  IP  PJ  4 xH  "3PJIAOX 
pj4oN.4^0tPp.rojDCVj^ 


p-  o' 
X z> 
O H- 


UJ  X X X O O ! 
XOOl-«XOCXO* 

h m a uj  or  ih  or  < 

Y UJ  UJ  IJ  r H o H 

I UJ  X X O CD  O c 

a in 

jh.t-hhJ-a.1 

D(AWWWWUJO 

OUJHHHHOffl 

xx<x<x<t<*»-i_i 

MOlIIIfflUJ 


: x r>  0 

■ I-  1/)K 
> O Z 4 


4 x-4 

PJ  CO 

fO  CVJ 


Z>  4 Z>  xO  Z)  Z>  f 

to  ^ CO  <D  to  ro  to  c 

ZIAZlAZIAZt 


I 4 UV  CO 
I • CVJ  .CVJ  • PJ  • Q 
CP  x-4  cO  H IA  UJ 

P.  4 p.  ro  h 

ro  ul 


H 


11 


CP  D N D 
I J)  1/J  H (/) 


CPZ>roz>^Z>UVO 

CPtoCPLOtPtor^o 
4 Z UV  Z UV  Z uv 


QC  _J  J rl  .CVJ 
< xH  * 

UJ  O P ^ V 

ac  z o fs,  r 

O < O xH 

U.  z ll.  11  II 


1 xO  • P~  • ao 
) xO  4 


• O'  • o • o *0. 

4 ^4  x-4  v£)  O P.  3T 

O uv  Q H (\J  H 


o cp 

HPjn4lAvOKCOO'Or4PJMUJP 
H H H r4  v ^ 


1 IL 

© 


Jt 


C_>oocvJt_>roc_>oOoorooiAOP»C_>4(_>o< 
UJP-UJP-UJCOUJ^UJxHUJUVUJxOU'cOUJcOUJ  4 i 
QCHn'HOCHtfPJQCPJQCHiyHQCHaf^QCHi 


© 


© 


69 


Figure  28  (cont'd) . Input  verification. 


SUMMARY  STATISTICS  FOR  COMBINATION  NO 


1 


< 2, 


3 


*) 


X = STATURE  Y = SHOULDER  HEIGHT  Z = CHEST  HEIGHT 

THE  MEAN  VALUE  = 1774.44  THE  MEAN  VALUE  = 1516.74  THE  MEAN  VALUE  = 

ST  * 0 OE  VI  AT  ION  = 66.47  ST*0  DEVIATION  x toO.55  ST'O  OEVIATION  = 


X X ( 

.695) Y, 

♦ < .372) Z 

4 

240.673 

SE 

= 24.542 

MULTIPLE  cork. 

COEF. 

Y = ( 

.312) X, 

♦ ( .709) Z 

* 

49.536 

SE 

= 16.442 

MULTIPLE  CORR. 

COEF. 

Z = ( 

. 164) X, 

♦ ( .69<,)Y 

♦ 

-54.023 

SE 

= 16.270 

MULTIPLE  CORR. 

COEF. 

SIMPLE  CORRELAT IONS-R( XY) 

= 

.925,  R< XZ) 

= 

.909,  R(YZ) 

= . 952 

N 

SUMMARY 

STATISTICS 

FOR  COMBINATION 

NO. 

2 

( 

3,  4,  5) 

X = SHOULDER  HEIGHT 

Y = 

CHEST  HEIGHT 

Z 

= WAIST  HEIGHT 

THE  MEAN  VALUE  = 

1516 . 7 A THE 

MEAN  VALUE  = 

1288.87 

THE  MEAN  VALUE  = 

ST'D  OEVIATION  = 

60.55  ST  *0 

DEVIATION  = 

54.72 

ST'D  OEVIATION  = 

X = ( 

.605) Y, 

♦ ( . 53  1)  Z 

4 

170.542 

SE 

= 16.022 

MULTIPLE  CORR. 

COEF. 

Y = ( 

.585) X, 

♦ ( . 3 5 7 ) Z 

4 

21. 277 

SE 

= 15.751 

MULTIPLE  CORR. 

COEF. 

Z = ( 

.483) X, 

♦ ( .336) Y 

4 

-98. 122 

SE 

= 15.273 

MULTIPLE  CORR. 

COEF. 

SIMPLE  CORRELATIONS-R(XY) 

s 

.952,  R (X Z) 

- 

.944,  R ( YZ) 

= . 934 
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SUMMARY 

STATISTICS 

FOR  COMBINATION 
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SE 
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SUMMARY 

STATISTICS 
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ST 

*0 

DEVIATION  = 

13.27 

ST  * O OcVIATION  = 

X = ( 

1.063) Y, 

♦ ( 

.291)2 

4 

-61.533 

SE 

x 14.266 

MULTIPLE  CORR. 

COEF  < 

Y = ( 

.264) X, 

4 ( 

.423 )Z 

4 

75. 015 

SE 

X 7.112 

MULTIPLE  CORR. 

COtF  i 

Z = < 

.121)  X, 

♦ ( 

.706)  Y 

4 

48.27 7 

SE 

= 9.189 

MULTIPLE  CORR. 

COEF, 

SIMPLE  CORRELATIONS-R(XY) 

= 

.763,  R ( XZ ) 

= 

.671,  R(YZ) 

= .775 

N 

SUMMARY  STATISTICS 

FOR  COMBINATION  NO. 

5 

< 

2,  12,  13) 

X = STATURE 

Y x 

HAND  LENGTH 

Z 

= FOOT  LENGTH 

THE  MEAN 

VALUE  = 

1774.44 

THE 

MEAN  VALUE  = 

193.18 

THE  MEAN  VALUE  = 

ST  * O OEVIATION  = 

66.  47 

ST 

*0 

DEVIATION  x 

8.85 

ST'D  OEVIATION  = 

X = ( 

3.064) Y, 

4 ( 

2. 27  9) Z 

4 

562.070 

SE 

x 43.950 

MULTIPLE  CORR. 

COEF, 

Y x ( 

. 043) X, 

4 ( 

. 428 ) Z 

4 

.180 

SE 

= 5.210 

MULTIPLE  CORR. 

COEF, 

Z = ( 

• 053)  X, 

4 ( 

• 7 0 9)  Y 

4 

41. 163 

SE 

x 6.7  03 

MULTIPLE  CORR. 

COEF 

SIMPLE  CORRELATIONS-R(XY) 

— 

.709,  R(XZ) 

- 

.705,  R(YZ) 

= .775 

N 

SUMMARY  ST 

ATISTICS 

FOR  COMBINATION  NO. 

6 

( 

1,  7,  8) 

X = WEIGHT 

Y = 

WAIST  BREADTH 

Z 

= WAIST  CIRC 

THE  MEAN 

VALUE  = 

179.79 

THE 

MEAN  VALUE  = 

322.05 

THE  MEAN  VALUE  = 

ST'O  OEVIATION  = 

20.70 

ST 

•0 

DEVIATION  x 

21.76 

ST'D  DEVIATION  = 

X = ( 

.471) Y, 

4 ( 

• 0 69) Z 

4 

-34.620 

SE 

= 14.950 

MULTIPLE  CORR. 

COEF 

Y = ( 

. 191) X, 

4 ( 

• 262 ) Z 

4 

50.105 

SE 

= 9.528 

MULTIPLE  CORR. 

COEF 

Z = ( 

.262) X, 

4 ( 

2 • 441 ) Y 

4 

74.275 

SE 

x 29.092 

MULTIPLE  CORR. 

COEF 

SIMPLE  CORRELATIONS-R(XY) 

= 

.685,  R(XZ) 

s 

.653,  R ( YZ) 

= . 888 

N 

Figure  29.  Summary  statistics  for  the  first  iteration. 
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= .958 
= .951 

= 150. 


907. 49 
63.96 
= .683 
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= .891 

- 150. 


70 


CO  00 


~r 


t 


• • • • 

yO  t4  t4  fO 

N CO  K CD 

CM  CM  O fO 


• • • • 

n.  in  ~t  <o 

r-  cm  o 

eg  M t4  cm 


X X 
O k— 
►H  J 
.d  Z 
X tjJ 


y»  ►-  in 

IlJ  O l-l 

I o «r 

O U.  X 

"0  ao 


O 

CX 

M 

O 


* * * * 

* * » » 

U?  CM  XS  CD 

• • • • 

in  CM  CM  N. 

f I 


CO 

0) 

3 

I— < 

cd 

> 

4J 

a 

td 

p 

p 

0) 

A 

cd 


in  cm 
~r  cd 
in  cm 


J-  H 

in  o> 


(X  z 

Id  id 
O —I 

3 O 
O Z 
X -*X 
(/)  X 

ro  cm 


-P 

•H 

CO  • 

-P  G 

a o 

d)  *H 
-P 
^ td 
3 P 
CO  CD 
-P 
M-l  -H 


-P 

CO 

■H 


-P 

CO 

P 

•H 


PI  IP 


go  uO 

O'  O' 


O 

CX 


o in 
in  a. 


o 

ex 


Ol  01 

K 

O CM  CM 


CO 

►- 

O O' 


O 

co 


in  -f 
■H  S. 
N.  N. 


<D 

p 

3 

•H 

Pu 


3 0 0 

in  z z 


71 


SUMMARY  STATISTICS  FOR  COMBINATION  NO 


1 


( 2,  3 


k) 


X = STATURE 

Y = SHOULOER  HEIGHT 

Z 

= CHEST  HEIGHT 

THE  MEAN 

VALUE  - 

1773.19 

THE 

mean  value  = 

1517.16 

THE  MEAN  VALUE  - 

1289.15 

ST  * D DEVIATION  = 

64.  91 

ST-0 

DEVIATION  - 

60.53 

ST-0  OEVIATION  = 

54.79 

X = ( 

.785) Y, 

♦ < 

.301)  Z 4 

194.619 

SE 

= 13.925 

MULTIPLE  CORR. 

COEF. 

= .977 

Y = ( 

.666) X, 

♦ ( 

.301) Z 4 

-52.604 

SE 

= 12.829 

MULTIPLE  CORR. 

COtF. 

= .9  77 

Z = ( 

. 387) X, 

+ < 

• 457 ) Y 4 

-91.109 

Sc. 

= 15.807 

MULTIPLE  CORR. 

COEF. 

= .957 

SIMPLE  CORRELATIONS 

-R(XY)  - 

.974,  R<XZ) 

s 

.951,  R(YZ) 

= . 952 

N 

= 149. 

SUMMARY  STATISTICS 

FOR  COMBINATION  NO. 

2 

( 

3,  4,  5) 

X = SHOULOER  HEIGHT 

Y = CHEST  HEIGHT 

Z 

= WAIST  HEIGHT 

the  mean 

VALUE  = 

1516.74 

THE 

MEAN  VALUE  = 

1 

288.87 

THE  MEAN  VALUE  = 

106  o. 39 

ST  * 0 OEVIATION  = 

60.55 

S T ® O 

OEVIATION  = 

54.72 

ST-D  DEVIATION  = 

49.45 

X = ( 

.605) Y, 

+ ( 

•531) Z 4 

170.542 

SE 

= 16.022 

MULTIPLE  CORR. 

COEF. 

= .964 

r = ( 

.585) X, 

♦ < 

.357) Z 4 

21.277 

SE 

= 15.751 

MULTIPLE  CORR. 

COtF. 

= .958 

Z ^ ( 

.483)  X, 

♦ < 

• 3 3 6)  Y 4 

-98.122 

SE 

= 15.273 

MULTIPLE  CORR. 

COEF. 

= .951 

SIMPLE  CORRELATIONS 

-R(XY)  = 

.952,  R(XZ) 

= 

.944,  R(YZ) 

= . 934 

N 

= 150. 

summary  statistics 

FOR  COMBINATION  NO. 

3 

( 

6 , 7 , 8 ) 

X = WAIST 

OEPTH 

Y = WAIST  9REAOTH 

Z 

= WAIST  CIRC 

the  mean 

VALUE  = 

235.41 

THE 

mean  value  = 

322.05 

THE  MEAN  VALUE  = 

907.49 

ST  * 0 OEVIATION  = 

20.  42 

ST-0 

DEVIATION  = 

21.76 

ST  * O OEVIATION  = 

63.96 

X = ( 

-.203) Y, 

♦ < 

• 3 41 ) Z 4 

-9.151 

SE 

= 9.577 

MULTIPLE  CORR. 

COtF. 

= . 3 83 

Y = ( 

-.211)  X, 

♦ < 

.361) Z 4 

43. 813 

SE 

= 9.772 

MULTIPLE  CORR. 

COEF. 

= .893 

Z = < 

1.531) X, 

♦ ( 1 

• 557 ) Y 4 

45.700 

Sc 

= 20.282 

MULTIPLE  CORR. 

COEF. 

= .948 

SIMPLE  CORRELATIONS 

-R(XY)  = 

.734,  R(XZ) 

at 

.878,  R(YZ) 

- .888 

N 

= 150. 

SUMMARY  STATISTICS 

FOR  COMBINATION  NO. 

4 

( 

9,  10,  11) 

X = BICEPS 

CIRC 

Y = ELBOW  CIRC 

Z 

= FOREARM  CIRC 

THc  MEAN 

VALUE  = 

316.26 

THE 

MEAN  VALUE  = 

278.09 

THE  MEAN  VALUE  = 

283.40 

ST-0  OEVIATION  = 

22.67 

ST-0 

DEVIATION  = 

13.1b 

ST-0  OEVIATION  = 

13.  52 

X = ( 

.482) Y, 

♦ ( 

. 876 ) Z 4 

-66.050 

SE 

= 13.915 

MULTIPLE  CORR. 

COtF. 

= .790 

Y = < 

. 055) X, 

♦ ( 

. 835 ) Z 4 

23.991 

SE 

= 4.701 

MULTIPLE  CORR. 

CUEF. 

= .934 

Z = ( 

. 099) X, 

+ ( 

• 82 7 ) Y 4 

22.184 

SE 

= 4.677 

MULTIPLE  CORR. 

COtF. 

= .938 

SIMPLE  CORRELATIONS 

-R(XY)  = 

.767,  R ( X Z ) 

— 

.783,  R(YZ) 

= .932 

N 

= 149. 

SUMMARY  STATISTICS 

FOR  COMBINATION  NO. 

5 

( 

2,  12,  13) 

X = STATURE 

Y = HANO  LENGTH 

Z 

= FOOT  LENGTH 

THE  MEAN 

VALUE  = 

1773. 19 

THE 

mean  value  = 

193.19 

THE  MEAN  VALUE  = 

272.19 

S T • 0 OEVIATION  = 

64.91 

ST'O 

DEVIATION  = 

8.88 

ST  * O OEVIATION  = 

11.  36 

X = ( 

3.220) Y, 

♦ ( 2 

•157) Z 4 

563.900 

SE 

= 41.325 

MULTIPLE  corr. 

COEF. 

= .771 

Y = ( 

•050) X, 

4 ( 

• 4 0 3)  Z 4 

-4.538 

SE 

= 5.134 

MULTIPLE  CORR. 

COEF. 

= .816 

Z = ( 

. 057) X, 

4 ( 

.688) Y 4 

38. 583 

SE 

= 6.707 

MULTIPLE  CORR. 

COtF. 

= .807 

SIMPLE  CORRELATIONS 

-R(XY)  = 

.734,  R(XZ) 

- 

.719,  R ( YZ ) 

= . 776 

N 

= 149. 

SUMMARY  ST 

ATISTICS 

FOR  COMBINATION  NO. 

6 

( 

If  7,  8) 

X - WEIGHT 

Y = WAIST  BREAOTH 

1 

: = WAIST  CIRC 

THE  MEAN 

VALUE  = 

183.15 

THE 

MEAN  VALUE  = 

321.56 

THE  MEAN  VALUE  = 

906. 04 

ST  * 0 OEVIATION  = 

20.28 

ST-0 

DEVIATION  = 

20.99 

ST-0  DEVIATION  -= 

61.66 

X = ( 

.528) Y, 

4 ( 

.092) Z 4 

-73.095 

SE 

= 12.056 

MULTIPLE  CORR. 

COEF. 

= .804 

Y = ( 

•305) X, 

4 ( 

. 22 3 ) Z 4 

64. 530 

SE 

= 9.161 

MULTIPLE  CORR. 

COEF. 

= .900 

Z = < 

. 521) X, 

4 ( 2 

• 1 8 3 ) Y 4 

110.032 

SE 

= 28.660 

MULTIPLc  CORR. 

COEF. 

= .385 

SIMPLE  CORRELATIONS 

-R(XY)  = 

.793,  R(XZ) 

= 

.761,  R ( YZ) 

= .079 

N 

= 149. 

Figure  31.  Recomputed  summary  statistics. 
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with  particular  attention  to  the  multiple  correlation  coefficients. 
The  best  possible  predictive  regression  equations  are  for  those 
combinations  having  a multiple  correlation  coefficient  which 
approaches  1.00.  With  such  a regression  the  estimated  or  pre- 
dicted value  will  be  essentially  identical  to  the  read-in  value 
unless  the  latter  is  in  error.  Conversely,  the  poorest  level 
of  prediction  is  for  those  regression  equations  where  the  multiple 
correlation  coefficient  approaches  0.00.  In  such  instances,  the 
predicted  value  for  a variable  will  be  the  group  mean  value  and 
all  values  over  CK  units  of  standard  deviation  from  the  mean 
will  be  flagged  as  aberrant  values. 

In  actual  fact  neither  extreme  of  the  multiple  correlation 
coefficient  normally  occurs.  Some  anthropometric  dimensions, 
such  as  long-bone  measurements,  have  relatively  high  inter- 
relationships (>.90)  whereas  measurements  of  the  head  and  face 
have  relatively  low  interrelationships  ( < . 2 0 ) . As  the  multiple 
correlation  coefficient  becomes  smaller,  the  estimated  value 
becomes  closer  to  the  mean  of  the  read-in  data  and  the  range  of 
acceptable  read-in  values  expands  with  CK  standard  deviations 
as  the  limit.  Conversely,  as  the  multiple  correlation  coeffi- 
cient becomes  larger,  the  estimated  value  is  less  dependent  on 
the  mean  of  the  read-in  data  and  the  range  of  acceptable  read-in 
values  contracts  (see  Table  1)  . 


TABLE  1 

RELATIONSHIP  OF  THE  MULTIPLE  CORRELATION 
COEFFICIENT  AND  ACCEPTABLE  DEVIATIONS 
FROM  THE  ESTIMATE 

R Aberrant  Value  Must  Exceed 


0.00 

+ 

3.5  SD 

0.30 

+ 

3.325 

0.60 

+ 

2.8 

0.70 

+ 

2 . 5 

0.80 

+ 

2.1 

0.90 

+ 

1.54 

0.95 

+ 

1.09 

1.00 

0.00 

Where  it  is  found  that  a variable  is  predicted  poorly  by 
the  other  two  variables  of  a combination  (that  is,  the  "R"  value 
is  low) , it  may  be  necessary  to  try  other  combinations  to  attempt 
to  improve  the  predictive  level  of  the  combination.  Even  when 
combinations  are  not  ideal,  however,  it  is  possible  to  use  them 
to  edit  the  data.  Poor  judgment  in  the  choice  of  combinations 
reduces  the  number  of  aberrant  values  called  out  but  does  not 
cause  acceptable  values  to  appear  aberrant. 

In  editing  data,  the  last  two  sections  of  the  printout  are 
of  principal  significance.  The  first  of  these  (Figure  32)  is 
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the  table  of  values  for  deviations  from  the  mean  for  each  variable 
by  half  units  of  standard  deviation.  These  values  range  from  -5 
standard  deviations  to  +5  standard  deviations,  an  extremely  broad 
interval  for  most  anthropometric  dimensions.  The  second  part  of 
the  editing  output  is  the  tabulation  of  the  aberrant  or  suspect 
values  by  subject  (Figure  33) . The  first  line  of  Figure  33 
contains  the  subject  number  followed  on  succeeding  lines  by  the 
combinations.  Each  line  contains  a separate  combination  with 
variable  number  and  name,  the  read-in  and  regression  values  as 
well  as  the  number  of  units  of  standard  error  of  estimate  these 
deviate  from  one  another  for  each  of  the  three  variables.  There 
is  such  a line  printed  out  for  each  subject  for  every  combination 
in  which  the  estimated  value  of  a variable  exceeds  the  read-in 
value  by  more  than  the  specified  CK  units  of  standard  error  of 
estimate.  Following  on  the  next  line  are  printed  out  the  vari- 
able names  called  for  in  the  list  associated  with  each  combina- 
tion. Below  each  variable  name  is  the  read-in  value  for  that 
variable  and  the  number  of  units  of  standard  deviation  it  lies 
above  or  below  the  variable  mean  value.  On  the  line  below  this, 
under  each  list  variable,  the  value  of  the  aberrant  or  suspect 
value  is  estimated  as  though  it  deviated  from  its  mean  as  many 
units  of  standard  deviation  as  the  list  variable  deviated  from 
its  mean  value. 

For  example,  as  can  be  seen  in  Figure  33,  subject  113  has 
been  flagged  as  having  a recorded  stature  of  1961  mm  (approximately 
6 Ft  3 in)  whereas  an  estimated  1711  mm  (approximately  5 ft  9 or 
10  in)  is  predicted  by  the  multiple  regression  equations  based 
on  measurements  of  his  shoulder  and  chest  heights.  The  difference 
between  his  recorded  and  estimated  stature  is  the  equivalent 
of  18  units  of  standard  error  of  estimate.  Scanning  the  pre- 
dicted values  of  stature  based  on  all  the  dimensions  in  the  list, 
it  is  obvious  that  the  recorded  value  is  much  too  high.  After 
studying  the  list  one  is  led  to  suspect  that  a transposition  of 
digits  has  occurred  and  that  the  recorded  value  of  stature 
should,  in  all  likelihood,  be  1691  and  not  1961.  As  a means 
of  double-checking  this  assumption,  we  note  that  for  this  subject 
shoulder  height  and  chest  height  (dimensions  which  correlate 
highly  with  stature)  have  recorded  values  which  are  approximately 
one  standard  deviation  below  their  means.  Turning  back  to 
Figure  32,  note  that  a stature  value  equivalent  to  one  standard 
deviation  below  the  mean  is  equal  to  1708  mm.  Thus,  our  assumption 
of  1691  would  appear  to  be  a reasonable  one.  A similar  process  is 
continued  for  each  of  the  aberrant  values  for  each  of  the  subjects 
tabulated  in  Figure  33. 

Editing  requires  common  sense  and  flexibility.  One  should 
always  remember  that  aberrant  values  may  be  either  errors  or  an 
abnormally  sized  or  shaped  individual.  Errors  are  generated  in 
three  ways:  a value  misread  in  measuring,  a misrecorded  value, 
or  a mispunched  value.  The  methods  of  making  the  measurements 
and  recording  the  data  should  be  reviewed  by  the  person  who  is 
going  to  do  the  editing.  The  editor  should  try  to  anticipate 
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those  errors  most  likely  to  occur.  Mispunched  values  are  easiest 
to  correct  since  there  is  no  judgment  to  be  made  if  you  have  the 
original  data  forms.  Sometimes  a value  is  recorded  so  that  it  can 
be  read  in  two  or  three  ways.  Often  one  value  will  be  very  close 
to  the  estimated  value  and  can  be  readily  corrected.  Values  may 
be  recorded  in  the  wrong  place;  two  or  more  may  be  reversed  or 
mixed  or  recorded  several  places  out  of  order  causinq  a number 
of  variable  values  to  be  one  or  more  positions  out  of  place.  The 
readily  identifiable  errors  are  those  which  occur  when  one  vari- 
able must  be  larger  than  the  other,  as  in  stature  and  shoulder 
height.  Where  it  is  suggested  that  two  values  may  be  reversed, 
an  estimate  of  one  based  on  the  other  in  any  combination  will 
be  of  no  use.  If  both  values  are  identified  as  aberrant,  the 
list  should  be  used  to  select  other  variables  which  correlate 
well  with  the  two  suspect  values. 

If  a whole  sequence  of  variables  are  flagged  as  aberrant, 
neither  the  estimates  derived  from  the  combinations  nor  the 
estimates  from  the  list  may  be  of  much  use.  In  such  instances, 
estimates  can  be  made  by  using  the  table  of  deviations  from  the 
mean  for  each  variable  (see  Figure  32) . Such  estimates  will 
be,  at  best,  gross  approximations.  Stature  and  weight  are 
useful  for  establishing  the  number  of  units  of  standard  deviation 
a particular  subject's  measurements  are  from  the  group  mean  values. 
These  variables  establish  the  overall  body  size  of  the  individual 
and,  in  general,  have  the  fewest  number  of  errors  in  a data 
set.  If  the  subject  is  more  or  less  of  average  height  but 
two  standard  deviations  above  the  mean  in  weight,  the  long-bone 
measurements  and  lengths  can  be  estimated  as  approximately  mean 
values  and  the  circumference  and  breadths  as  the  mean  plus  1.5  - 
2.5  standard  deviations  above  the  mean.  The  values  so  derived 
will  often  provide  an  adequate  clue  to  sort  out  sequential 
data  which  have  been  misrecorded.  A list  of  the  approximate 
variable  values  can  be  moved  up  and  down  beside  the  data  blank 
until  a good  match  is  found. 

Sometimes  a number  will  be  recorded  incorrectly.  One  or 
more  digits  may  be  wrong  or  out  of  order.  If  a value  is  estimated 
as  830  and  read  in  as  735,  a good  guess  is  that  835  is  the  actual 
value.  If  a value  is  estimated  at  585  and  read  in  as  509,  it  is 

a good  guess  that  590  is  the  actual  value.  Values  misread  in 

measuring  are  of  two  types:  the  measurements  were  not  made 
correctly  or  the  measuring  instrument  was  not  read  correctly. 

This  program  does  not  detect  mismeasurements  well,  but  if  for  a 
particular  variable  suspect  values  occur  for  several  subjects 
measured  at  the  same  time  or  by  the  same  individual,  it  is  a 
good  idea  to  check  for  an  equipment  problem  or  a change  in  method 
of  measurement.  An  instrument  which  is  not  marked  in  hundredths 
is  typically  the  cause  of  some  measuring  errors  which  creep  into 
the  data.  For  example,  if  most  values  are  between  800  to  899 

with  only  a few  above  and  a few  below,  then  values  beginning  with 

7 or  9 , when  suspect,  may  be  corrected  to  begin  with  8. 
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While  the  EDIT  program  can  sort  through  a massive  amount  of 
data  to  locate  and  identify  deviant  values,  it  cannot  make 
judgments  about  those  aberrations.  It  requires  a knowledgeable 
professional  to  decide  whether  an  aberrant  value  is  in  fact  an 
error,  to  track  down  its  source,  if  practicable,  and  to  decide 
ultimately  whether  the  value  can  be  corrected  or  should  be 
deleted . 

The  EDITing  program  is  independently  valuable  but  its 
optimum  usefulness  is  achieved  in  conjunction  with  the  previously 
described  XVAL  program  which  serves  to  identify  the  grosser 
errors,  leaving  to  EDIT  the  task  of  fine  sorting. 
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APPENDIX  A 


COMPUTER  PRINTOUT  OF  THE  XVAL  PROGRAM 
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PROGRAM  XVAL77 

*< INPUT, OUT  PUT, TAP£5  = INPJT , TAPt&  = OUT  PUT, TAPE9, PUNCH, TAPE2, TAPE  3) 
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FILE  NAMES  MODE 

TAPE5  MIXED  READS  63  91  95  99  101  106 

TAPE6  MIXED  WRITES  73  81  92  102  109 
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STATISTICS 

PROGRAM  LENGTH  371B  249 

CM  LABELE0  COMMON  LENGTH  5415B  2929 
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PROGRAM  LENGTH  7B  7 
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COMPUTER  PRINTOUT  OF  THE  EDIT  PROGRAM 
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INLINE  FUNCTIONS  T*PE  ARCS  DEF  LINE  REFERENCES 

A HA  XI  REAL  0 INTRIN  173  174  175  197  198  199 

AMIN1  REAL  0 INTRIN  194  195  196 

MAXO  INTEGER  0 INTRIN  72 

SQROOT  REAL  1 SF  49  173  174  175  197  198  199 
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STATISTICS 

PROGRAM  LENGTH  13478  743 

CM  LA8ELE0  COMMON  LENGTH  134048  5892 
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STATISTICS 

PROGRAM  LENGTH  78  7 

CM  LA8ELE0  COMMON  LENGTH  31348  1628 


